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1.  INTRODUCTION  AND  SUMMARY 

1.1  INTRODUCTION 

The  scattering  of  an  electromagnetic  wave  by  a dielectric  wedge 

is  an  exceptionally  difficult  problem  which  at  the  present  time  has  no 

(1)  (2) 

known  analytic  solution,  though  Radlow,  in  the  special  case  of  a 
right-angle  wedge,  was  able  to  use  a double-Laplace  transform  method  that 
had  been  found  successful  in  treating  the  quater-plane  problem  to  analyze 
the  wedge  exactly.*  All  other  treatments  yield  approximations  of  one  kind 
or  another.  An  approximate  numerical  solution  for  the  diffraction 
coefficient  has  been  given, ^ but  the  shortcoming  of  this  sort  of  approach 
is  that,  although  engineering  and  other  applications  eventually  require 
numerical  calculations,  the  wedge  field  is  a compound  of  several  different 
components,  such  as  specularly  reflected  and  refracted  rays,  surface  waves, 
lateral  waves  and  tip  diffraction,  and  a numerical  calculation  of  the  com- 
posite effect  throws  very  little  light  on  the  real  nature  of  the  response. 

An  analytic  solution,  even  though  not  exact,  is  highly  desirable  in  order 
to  enable  the  different  components  of  the  solution  to  be  separated, 
recognized  and  understood.  A study  of  the  functional  form  of  the  result 
can  yield  considerable  insight  into  the  rather  complex  properties  of  the 
wedge,  and  In  particular  the  behavior  near  shadow  boundaries,  as  well  as 

the  effect  of  varying  such  parameters  as  wedge  angle  or  dielectric  constant. 
(4) 

Balling,  in  his  Licentiate  thesis,  provides  an  analytic  treatment 
based  on  tracing  the  multiply-reflected  rays  inside  the  wedge.  In  two  later 
papers, he  applies  his  results  to  evaluate  the  role  of  surface  fields 

Radlow's  results  are  too  involved  for  use,  and  Kuo  and  rionus^11^ 
and  Kraut  and  Lehman'1"'  have  attempted  alternative  formulations.  The  latter 
authors  in  fact  find  Radlow's  analysis  to  be  erroneous.  In  any  case  the 
method  seems  incapable  of  generalization  to  other  angles. 
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and  lateral  waves  in  the  radiation  from  a wedge  excited  by  a line  source 
inside  the  wedge.  All  these  treatments  utilize  formulas  for  successive 
ray  reflections,  and  Borovikov^  ^ has  also  obtained  a recurrence  formula 
for  the  calculation  of  the  terms  for  an  edge  wave,  provided  none  of  the 
edges  are  located  near  shadow  boundaries.  A somewhat  different  approach 
by  Tricoles  and  Rope^  based  on  the  penetration  of  the  field  in  a dielec- 
tric slab  calculates  the  approximate  effect  of  a hollow  dielectric  wedge 
made  by  two  slabs  meeting  at  a bevelled  edge.  The  application  is  to 

radomes,  where  the  hollow  structure  is  used. 

(9) 

Kaminetzky  and  Keller,  in  the  most  recent  published  results, 
analyze  the  wedge  in  two  cases  in  which  the  tip  diffraction  can  be  cal- 
culated approximately  as  the  first  term  in  an  expansion  of  the  field  in 
terms  of  a small  expansion  parameter.  The  two  cases  they  treat  are 
i)  a small  difference  in  the  dielectric  constants  of  the  materials  inside 
and  outside  the  wedge,  and  ii)  small  wedge  angle.  The  latter  complements 
an  earlier  treatment*'11^  where  the  wedge  angle  is  nearly  180°,  and  the 
expansion  is  relative  to  the  reflection  of  a wave  at  a plane  dielectric 
surface . 

The  wedge  problem  is  an  important  one  in  at  least  two  areas.  The 
first  concerns  radar  reflections  and  BMP  pulse  response  from  dielectric 
objects  which  may  be  in  free-space  or  else  buried.  The  second  concerns 
the  use  of  the  geometric  theory  of  diffraction  (GTD)  to  calculate  the 
radiation  properties  of  antennas  and  other  reflectors.  For  example,  a 
dielectric  support  may  be  used  to  locate  a subreflector  or  splash-plate 
and  the  effect  of  the  dielectric  in  modifying  the  diffracted  fields  from 
metallic  edges  needs  to  be  taken  into  account.  The  radiation  from  most 
composite  reflectors  can  be  represented  in  terms  of  specular  reflections 
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and  edge  effects,  and  the  advent  of  a new  class  of  formulas  involving 
dielectric  wedges,  either  on  their  own  or  in  contact  with  a metallic  surface, 
can  be  expected  to  extend  the  range  of  application  of  the  GTP  method. 

The  shadow -boundary  terms  are  an  important  part  of  the  toal  solution 
and  should  be  shown  explicitly. 

Although  the  geometry  of  the  arrangement  is  quite  simple,  nevertheless 
the  existence  of  propagating  waves  inside  and  outside  the  wedge,  with 
different  velocities,  makes  the  problem  a very  formidable  one. 

1 • - SUMMARY 

Three  separate  investigations  were  made.  The  first,  based  on  a 
spectral  decomposition  method,  was  the  one  on  which  the  contract  proposal 
was  originally  based.  It  turned  out  to  have  a built-in  flaw,  which  was 
not  discovered,  however,  until  rather  late  in  the  analysis.  As  a con- 
sequence the  method  was  abandoned.  Just  at  that  time,  attention  was  drawn 
to  two  Soviet  papers,  both  of  which  claimed  to  have  solved  the  wedge 

problem  rigorously  and  completely,  though  by  completely  different  methods. 

< i jf 

The  first  paper,  by  Zavadskii,  was  examined  and  found  to  be  faulty. 
Efforts  made  to  correct  the  errors  were  to  no  avail.  However,  a semi- 
trivial  result  emerged  from  these  studies  which,  in  retrospect,  could  have 
been  found  by  elementary  methods.  The  second  paper,  by  Aleksandrova  and 
Khicnyak  , also  turned  out  to  be  faulty,  though  the  complex  analysis 
involved  a great  deal  of  study  to  pin-point  exactly  where  the  errors 
occurred.  It  was  not  found  possible  to  solve  the  problem  by  their  method 
with  the  errors  corrected. 

It  appears  that  at  the  present  time  there  is  still  no  valid  and 
rigorous  solution  to  the  wedge  problem,  though  several  approximations,'  ' 
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particularly  Balling’s,  have  appeared  in  the  literature,  and  should 
provide  useful  results  in  practical  applications. 
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DIFFRACTION'  BY  A DIELECTRIC  WEDGE:  SPECTRAL  ANALYSIS 


ABSTRACT 

A spectral  analysis  of  the  field  inside  and  outside  a dielectric  wedge 
is  made,  and  field  matching  in  the  presence  of  an  axial  incident  wave  is 
made,  both  on-axis  and  on  the  wedge  faces.  The  method  involves  unknown 
spectral  functions  among  which  a relation  (14)  is  imposed  in  order  to 
lead  to  Wiener-Hopf- 1 ike  equations  (21)  and  (22).  However  an  examination 
of  the  pole  locations  of  these  equations  shows  that  they  self-generate 
an  unending  pole-sequence,  which  makes  the  method  unusable.  No  modification 
of  the  method  that  would  render  it  valid  appears  possible. 
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Pespite  its  simple  geometry,  the  problem  of  a dielectric  wedge  excited 
I"  an  incident  wave  is  an  extremely  difficult  one,  and  no  exact  solutions 
are  known.  Various  features  l\a\e  been  treated  by  different  researchers, 
including  an  examination  of  the  field  near  the  tip.  specular  internal  and 
external  reflection,  and  properties  when  the  dielectric  constants  of  the 
wedge  and  t he  surrounding  medium  are  nearly  the  same. 

In  this  paper  we  attempt  an  exact  formulation  from  which  useful  results 
can  be  extracted.  In  view  of  the  extreme  difficulties  encountered,  the 
total  wedge  angle  has  been  restricted  to  less  than  90°  (other  features 
apparently  enter  when  the  wedge  angle  is  larger  than  this),  and  the  analysis 
is  further  confined  to  symmetrical  incidence  of  a plane  wave  with  parallel 
po lar i ;a t ion . 

i'oumiiluion 

Initial  attempts  to  provide  a formulation  in  terms  of  cylindrical 
functions  were  discouraging,  and  eventual  l'  an  expansion  of  the  fields  as 
an  integral  of  plane  waves  was  selected  as  more  promising.  However,  a number 
of  features  had  to  be  closely  constrained  in  order  to  permit  field  matching 
in  a constructive  way  and  this  accounts  in  part  for  t he  form  taken  by  the 
subsequent  anal' sis. 

figure  J.l  shows  the  eonf igurat ion  in  which  a dielectric  wedge,  of 
angle  -'o  and  refractive  index  n is  immersed  in  a medium  of  refractive  index 
n^.  The  latter  would  normally  correspond  to  free  space,  but  it  turns  out 
to  be  essential  to  allow  both  media  to  possess  a small  loss  term.  This  will 
be  taken  to  the  zero  in  the  limit,  but  for  much  of  the  analysis  it  has  to  be 
retained,  since  it  determines  the  a 1 1 - important  feature  of  whether  certain 
singularities  are  inside  or  outside  integration  contours. 


l-'ig.  2.1  Geornet ry  o 1'  p io  1 oc trie  Wedge 

1-igurc  2.1  shows  a rectangular  coordinate  system  with  the  x-axis 
bisecting  the  wedge  and  the  z-axis  along  its  tip.  A cylindrical  coordinate 
system  p,  41,  2 is  also  erected,  with  4>  = ±a  at  the  interfaces. 

The  total  region  is  divided  into  four,  as  follows: 

Region  I 01  < < tt 

II  0 < <(>  < a 

I ' it  < <j>  < 2tt  - cx 


In  view  of  the  chosen  symmetry  of  the  incident  wave  the  fields  in  Regions  1 
and  11'  will  be  mirror  images  of  those  in  1 and  11.  and  can  be  allowed  for 
by  replacing  (.it  - 40  by  |tt  - 4>|  and  4>  by  j,"J  respectively.  However,  this 
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does  not  automatically  ensure  field  matching  on  the  axis,  since  II 

P 

proportional  to  - (.  1 / p ) L. II  / will  in  general  be  discontinuous  on  axis  due 
to  the  presence  of  the  modulus  terms  in  the  azimuth  angle.  However,  by 
requiring  3l:.,/34>  to  be  zero  on  axis,  corresponding  to  II.  = 0 there,  as 
required  by  symmetry,  the  necessary  continuity  is  obtained.  1'he  reason 
for  this  choice,  rather  than  for  the  more  usual  one  of  selecting  functions 
of  cos  41  • which  are  automatically  symmetrical  on  axis,  stems  from  the 
form  taken  by  the  fields  at  41  = -o.  . To  enable  matching  in  a form  that  is 
mathematically  tractable  it  seems  necessary  to  work  with  the  above 
formulat ion . 

-•5  l-ir.l.H  l'XPUI'SS I ONS  IN  Kltl  1 ON  I 

Throughout  this  paper  we  take  w as  a spectral  variable  and  build 

solutions  in  region  l from  the  basic  plane  wave 

*»  *»  1 

-p[  (w“-n“)  2sin(.4>-cx)  - jwcos  (4>-ctl  ] (.11 

E.  = c ° 

In  this  expression  p is  the  radial  distance  from  the  tip.  normalized  with 

respect  to  the  free-space  wave  number  k . i.e.  p ;;  k p where  p is  the 

O h O O o 

actual  radial  distance.  1’he  exponential  is  put  in  this  form,  one  of  many 

possible,  in  anticipation  of  eventual  field  matching  at  41  = a . 

As  mentioned  earlier,  n contains  a small  loss  term,  and  we  put 

0 

no  = ~ ln0l  IM«)  (.2) 


Equation  (.1)  has  branch  cuts  at  w = ±n0»  and  in  order  to  avoid  divergent 

o ->  9 

waves  at  infinity  that  branch  of  (,w“  -n-)  has  to  be  chosen  to  give  damped 
waves  as  w * This  leads  as  one  possibility  to  the  branch  cuts  shown  in 

Figure  J. 2a,  but  there  is  some  latitude  as  to  where  the  cuts  go. 


w- PLANE 


•Lg.  2.2  Branch  Cuts  in  w-plane 


Since  a < $ < tt  in  region  I it  is  found  that  if  w = Re*^  as  R 03  that 


0 < 3 < a is  a possible  range,  whilst  if  w = Re*  as  R -*■  -w,  then 


-a  < 3 < 0 is  a possible  range  for  3 . It  turns  out  from  other  consider- 


ations that  a branch  cut  at  3 = -a  is  needed  at  w and  this  leads  to 


the  modified  branch  cut  shown  in  Figure  2.2b.  For  the  time  being  we  take 


w exactly  on  the  real  axis  and  get,  for  the  scattered  electric  field  i 


region  I 


2 2 L 

E.  = [ f (w)e-pt  (w“~nG)  "sin  (<J>-a) -jwcos  (<(>-u)  ] 


dw  (3) 


with  f(w)  an  undetermined  spectral  function.  Before  proceeding  we  need  to 
check  both  for  completeness  and  performance  at  infinity.  As  far  as  the 
latter  is  concerned,  this  is  covered  by  the  convergence  of  (3)  as  p 
due  to  the  choice  of  branch  cut.  However,  individual  plane  wave  components 
need  not,  in  isolation,  satisfy  the  radiation  condition,  which  is  a restraint 
on  the  total  field  only.  Specifically,  when  — | n | < w < |n  | we  can  put 
w = In  | cos  0 with  0 < 0 < tt,  and  an  individual  term  in  (3''  becomes  like 
a plane  wave*  travelling  at  angle  (J)  = it  + a - 0.  When  0 < 0 < a this  gives 
waves  at  angles  greater  than  tt,  which  would  not,  in  the  absence  of 
attenuation,  be  acceptable  in  isolation.  The  fact  that  (3)  is  integrated 
to  +<»  ensures  in  any  case  the  correct  performance  of  the  total  field,  as 
is  eventually  confirmed  when  the  fields  are  found. 

As  far  as  completeness  is  concerned,  it  might  appear  that  a term 

9 

similar  to  (3)  but  with  the  sign  of  (w“  -n~)~  reversed  would  also  be 

o 

needed,  since  such  a term  would  also  satisfy  the  wave  equation.  However, 
for  w"  > |n  | “ the  wave  would  diverge,  whilst  for  -In  I < w < In  I the 
wave  components  would  represent  incoming  waves  at  angles  from  it  + a to 
— tt  *•  a;  not  only  are  they  physically  irrelevant  when  taken  in  isolation, 
but  they  turn  out  to  have  i nwa rd  attenuation  when  w is  real,  i.e.  they 
arc  incoming  waves  from  infinity,  and  form  no  part  of  the  solution. 


An  incident  field  on  axis  of  unit  amplitude  is  simply 
whence  we  get  the  total  fields  in  a < <p  < it  , 


-lpn  cos  >t> 
o 


The  IS  term  in  n ensures  that  these  waves  arc  attenuated  outward  when  w 
is  real.  ° 


Ez  = o"jpnoCOS'  + j e"Pt(w  -no)?sin(<J!-a)-jwcos(^,x)lf(w)dw 

— .v> 

(4) 


- = -jn  siml'  e‘*pnocos* 

p 4<j>  o 


i o V 


+ j e_pl(w  'no>  ^in(^)-jwcos(^-u)]f,w)  [lw2.n2)  l'Cos(^-« 


+ jwsin  ( vfi-a ) ] dw 


(5) 


(The  latter  expression 


is  proportional  to  U 


2.4  FIELD  EXPRESSIONS  IN  REGION  II 

The  analysis  follows,  in  part,  the  considerat ions 

2 ? *>  2 2 \ 

outlined  in  section  2. 3 with  (w  -n  ) replacing  (w  -n  ) ; but 

o 

there  are  some  important  differences.  In  region  II  the  waves 

2 2 \ 

corresponding  to  reversing  the  sign  of  (w  -n“)  are  not  only 
acceptable,  from  the  point  ol  view  of  convergence  at  infinity; 
they  are  also  required,  and  in  fact  represent  those  waves 
incident  in  region  II  from  region  IT',  where  they  are  generated 
by  refraction  at  the  boundary  at  i)>  = -a  and  cross  the  axis  into 
region  II.  Accordingly  we  get,  for  0 < 0 < a 


l\> 

Ez  = | e~pl(w  -n2)  Jsin(a-4.)-jwcos(a-4')  ]g^  (w)dw 


f CP  ( (w2-n 2 ) 'Jsin(a-i(i)  + jwcos  (a-<|> ) 1 g (w)  dw  (6) 


. 3E 

1.  7. 
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f e-pt(w2-n2)SSin(«-*) -jwcos  («-J.)]gi(w)  [ . (w2_n2 , Scos  (a_* , 


— 1\> 


- jwsin  ) ldw  + 


integrals  represent  equal  downward  and  upward  waves  respectively 
within  the  wedqe.  The  contour  tor  g^  is  the  real-axis,  but  that 
for  g is  not,  since  its  integrand  does  not  converge  (for  0 if  • a) 

at  w +•"  on  axis.  In  fact  we  need  w - Re  ,lX  and  -Re  ,a  as 

R > in  order  to  give  an  integrand  which  has,  for  the  range 
O*'  if  a,  a non-positive  real  part  t o the  exponent.  We  shall 
later  choose  this  contour  to  correspond  t o a purely  imaginary 

exponent  at  if  a,  but  it  is  sufficient  .it  this  stage  to  note  that 

suitable  contours  exist , and  that  the  left-hand  branch  cut  must 
be  t i 1 t mi  upwards  (.it  Infinity)  by  an  angle  a to  accommodate  this 
fea t ure . 

I'lllt,  0 MATCH  INC  ON  AXIS  AT  ,f  n 

1 

As  mentioned  earlier,  the  field  in  1’  is  obtained  from  I 
by  replacing  (11-if)  by  | 11  - if  | • At  if  n this  does  not  in  any 
way  attect  , which  is  therefore  continuous  at  the  boundary;  but 
i .'if  changes  sign  when  operating  on  | n - if  | as  if  goes 

Since  is  proportional  to  11  , which  in  the 

J 


through  n . 
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symmetrical  case  vanishes  at  4>  = tt  anyway,  it  is  clear  that 
matching  on  the  axis  requires  dE^/2<\>  = O at  if  = it.  From  (5) 
with  4>  = it  we  thus  get 


0 = 


e ^ Uf  (w)  [ (w~-n2  ) 'cosa-  jwsina]  dw  0 < p < 00  (8) 

o 


where 


U 


, 2 2,  Vj 

(w  -n  ) s 


ina+ jwcosa 


(9) 


Now  3U/3w  = j (w2-n2)  J [ (w2-n^) ^cosa- jwsina] , and  if  we  define  a 
new  function  F by 


f (w) (w2-n2 ) * = F(U) 


(10) 


then  (8)  can  be  written 
w=+>» 

| e'pUF(U)  dU  = 0 0 < p < 00  (11) 

w=-°° 

It  is  clear,  therefore,  that  (11)  can  be  satisfied  by  choosing 
F (U) , qua  function  of  U,  to  be  free  of  singularities  in  the 
region  Re  U > 0.  Moreover,  convergence  of  (11)  as  p -*■  0 requires 
F (U)  = 0 ( | U | 1)  as  U approaches  infinity  along  the  contour 

in  the  U-plane  corresponding  to  the  real  axis  of  w.  This  contour 
is  shown  in  figure  2.3;it  starts  at  U = °°e  comes  in  close 

to  the  imaginary  U-axis,  rises  to  a maximum  at  U = nQ  and  then 
drops  to  the  right  before  rising  asymptotically  to 
U = «>e^n//2  . The  closeness  of  approach  to  the  imaginary  axis 

is  proportional  to  6.  Also  shown  on  the  contour  is  the  line 
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Fig.  2.3  Contour  C^  in  the  U-plane 

U = jn  cos(O-a)  corresponding  to  w =n  cos  0.  This  intersects 
° o 

the  contour  at  0 = and  (almost)  at  0 = a;  and  between  these 
positions  the  line  w = nQCos  0 corresponds  to  values  of  U to 
the  right  of  the  contour.  Since  we  shall  later  be  concerned 
with  poles  at  various  positions  along  w =n^cos  0 it  is  clear 
that  poles  for  0 v 0 < it/2  cannot  be  permitted.  It  is  also  clear 
that  F (U)  free  of  singularities  for  Re  U < 0 is  unnecessarily 
restrictive;  it  is  only  necessary  for  it  to  be  analytic  between 
the  contour  C.  and  infinity. 
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It  should  be  noted  that  F(U)  analytic  qua  U is  not  the 

same  as  f (w)  analytic  qua  w.  In  fact  such  a function  as 

F (U)  = U is  clearly  analytic  in  U but  because  of  (9)  it  has 

branch  cuts  at  w = +n  . 

— o 

Around  the  point  w = n^cosln+iH  we  have  U = jnQcosii)  which 

is  therefore  symmetric  about  tp  = O.  Accordingly,  if  the  left 

hand  part  of  the  real-axis  contour  in  the  w-plane  is  swung  around 

the  point  w = n^cosa,  in  fact  to  be  asymptotic  to  the  line 
“ i 2 ot 

w = r >ve  , the  corresponding  contour  in  the  U-plane  can  be  made 
to  collapse  to  a curve  doubling  back  on  itself.  Figures  2. 4a  & 2.  4b 


Fig.  2.4  Contours  in  the  w and  U-planes 
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outline  the  corresponding  contours  and  show  that  there  is  a 
region  in  Im  w ^ 0 which  is  outside  the  region  which  gives  the 
contour  C^.  Since  the  contour  of  Figure  2.4b  lies  within  the 
contour  it  is  apparent  that  (8)  is  satisfied  for  F(U) 
analytic,  qua  U,  between  and  infinity. 

2.6  FIELD  MATCHING  ON  AXIS  AT  | = 0 


In  the  same  way  as  described  in  section  2.5  we  require 

3E  /34>  to  be  zero  at  4>  = 0.  If  we  write 
z 


W,  W* 


2 2 \ 

+ (w  -n  ) sina- jwcosot 


(12) 


then  (7)  gives 


W = +"’ 


0 = 


, w 2 2,S  -pW ...  f . w 2 -pW* 

g^(w) (w  -n  ) e dW  - g2 (w) (w  -n  ) e dW* 

w=— w w=C 

(13) 

0 < p < °° 


The  contour  w real  gives  the  contour  C2  for  W,  as  shown  in  figure 
2.  5.  The  contour  C for  w in  the  g0  integral  is  somewhat  at 
choice.  The  larger  Im  w is  on  this  contour  the  larger  is  the 
real  part  of  W*.  There  is  a lowest  contour,  shown  in 
figure  2.6  for  which  Re  W*  = 0 and  C cannot  be  taken  below  this 
if  the  second  integral  in  (13)  is  to  be  convergent.  A little 
later  it  will  be  seen  that  we  wish  to  make  the  two  contours  in  the 
W and  W*  plane  as  close  as  possible,  so  we  choose  w = C to  give 
W*  imaginary  on  it.  If  we  further  choose  a new  function  G such  that* 


* 

This  choice  is  not  mandatory,  but  the  method  of  this  section 
is  dependent  on  it. 


w 


then 


Fig.  2.5  Contour  C,  in  the  W-plane 


g^(w)(w~-n~)!i  = G(Vtf),  g^lwliw^-n  ) ' - G (W* ) 


(13)  becomes 


f G (W) e pWdW  +• 


G(W*)e”pW  dW*  =0  O < p < -v'  ( 


Now  in  (15)  both  W and  W*  are  functioning  merely  as  dummy 
variables,  so  in  the  second  integral  we  can  replace  W*  by  W 


gj  - PLANE 


Fiq.2.6  Contour  C in  the  w-planc 


Moreover  the  contour  C,  can  be  joined  to  the  imaginary  axis  by  ar 
at  infinity  (which  give  zero  on  integration),  to  produce  the 
closed  contour  C^  of  f igure  .1 .7  cons  ist ing  of  (1)  the  imaginary 
axis,  (2)  arcs  at  +_  j>" , (3)  the  contour  C-  of  figure  2.5. 

Equation  (15)  can  now  be  written  simply 


* rW 

G (W)  e v dW  0 


0 v a v iv 


and  it  was  in  order  to  achieve  this  form  that  g^  and  u.  were 
chosen  to  be  related  as  in  (14).  Equation  (It')  is  satisfied  if 
G(W)  qua  W is  analytic  in  C..  Also  shown  in  figure  '.7is  the  line 
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2.7  FIFLD  MATCHING  AT  > = a 


„ , . -jn  pcosa 

From  (4)  the  incident  field  at  i)>  = a varies  as  e o 

and  before  we  can  field  match  it  is  necessary  to  put  this  in  the 

form  of  an  integral  compatible  with  the  others.  This  is  achieved 

through  writing 


e~ jn0pcosa 


CO 


ejpW 

w+n  cosct 
o 


dw 


(17) 


a result  that  is  readily  verified  by  deforming  the  contour  into 
Im  w > 0,  whence  the  pole  at  w = -nQCosa  provides  the  necessary 
residue.  This  is  only  valid  when  a < n/2,  but  since  we  are  later 
restricting  a to  be  not  greater  than  tt/4  the  constraint  is  of 
no  consequence.  Accordingly,  with  4>  = a , equating  (4)  to  (6) 
gives 


2 rr  j 


Jpw 


w+n  cosa 
o 


dw 


CO 

+ f(w)e-*pw  dw 

— CO 


g.  (w)  e^ pW  dw 


+jg2(w)e^pWdw  (18) 

C 

Before  proceeding  further  it  is  necessary  to  deform  C to  the 
real  w-axis.  This  could  not  be  done  before  because  with  the 
more  general  exponent  for  0 < < n the  integrand  was  not 

convergent;  but  with  e-^W  in  the  integrand  it  can  be  done, 
together  with  non-contributing  arcs  at  +00.  Now  if  g^  (w)  has 
any  singularities  between  the  real  w-axis  and  the  contour  C 
they  will  give  rise  to  a contribution  on  collapsing  C to  the  real 


J - L 6 


axis.  To  prevent  this,  write  y.,  = g,  + g where  g^  is  that 

part  of  y7  carrying  singularities  (poles  or  branch  cuts)  between 

C and  w real,  and  is  free  from  singularities  above  C.  Then  by 

deforming  C to  +j-v  we  see  that  f g (w)e^tWdw  = 0 and  (18)  can 

J s 

C 

be  written 


e3  Pw  [ f (w)  _g  (w) -g  (w)  + ( 1/2tt  j ) (w+n  cosa)  1 ] dw  = 0 


i 


0 < p 


(19) 


If  we  define  a plus  function  h (w)  to  be  free  of  sing 
for  im  w > 0 then  (19)  can  be  written  in  the  form 


f (w)  + TT— r [ 


2 n j w+n  cosa 
o 


- h (w) ] = g. (w)  + g2 (w) 


(20) 


Since  w = -n  is  above  the  C-contour  the  presence  of  branch  cuts 
at  this  position  doesn't  affect  g (w)  . On  putting  f and  g in 
terms  of  F and  G through  (10)  and  (14) , equation  (20)  becomes 


F (U) 


/ 2 2.  S 

(w  -n  ) 2 
o 


2 7T  j 


w+n  cosa 
o 


, , , , G (W) +G (W* ) 

- h (w)  ] = 5 

(w“-n“ ) J 


(21) 


Here,  G(W*)  comes  from  G(W*)  by  dropping  those  singularities, 
if  any,  that  occur  in  the  w-plane  between  C and  the  real  axis. 
If  a similar  process  is  used  to  match  at  $ = a the  following 
relationship  is  obtained 


n sina 

-FlU>  * ~2”n  1 


w+n  cos*  - ~ G'M>  - G<W*»  <22> 


o 
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where  k+  is  a similar  function  to  h + . Equations  (21)  and  (22) 
are  the  basic  relations  from  which  G and  F can,  in  principle, 
be  found,  and  which  then  determine  the  fields  via  (4)  and  (6). 

2.8  FOLE  CONTRIBUTIONS  TO  THE  SOLUTION 


Equations  (21)  and  (22)  both  exhibit  a pole  at  w = -nQcosa. 
This  is  at  a position  where,  if  g9 (w)  (or  G(W*))  possessed  a 
pole  it  would  have  to  be  subtracted  out  to  form  (w) . Hence 
G(W*)  in  (21)  and  (22)  cannot  contribute  such  a pole,  and  only 
F ( U)  and  G(W)  can  do  so.  That  both  must  do  so  follows  at  once 

from  the  different  ways  in  which  F,  G and  the  pole  enter  into 
the  two  equations. 

Let  U and  W be  the  values  of  U and  W at  w = -n  cosa. 

u U Q 

It  is  easily  found  that 

Uc  = -jnQCOs2a  (23) 

W = jncos(0  - a)  where  ncos0  = n cosa  (24) 

o Or  O O 


Since,  from  (24),  a < ©o  < tt/2  it  can  be  seen  from  figures  2. 3 and 
2.7  that  Uq  and  Wq  are  in  regions  where  F and  G may  possess 
poles,  and  so  we  put,  in  the  neighbourhood  of  these  poles 


A 

U-U 

o 


B 

W-V* 

o 


F(U) 


G ( W) 


(25) 


r 


2-18 


where  A and  B are  constants  to  be  determined.  The  residues  at 
w =-nQcosa  are  readily  found  from  (25)  and  the  definitions  of 


U and  W in  terms  of  w,  and  give 


_A „ A/2  jcosa 


U-U 


u-u 


w+n  cosa 
o 


B/2jcosa 

w+n  cosa 
o 


2nsin9 


n sina-nsin0 
o o 


(26) 


(27) 


Hence,  to  match  the  poles  at  w = -n^cosa  in  (21)  and  (22), 
since  they  cannot  be  provided  in  that  position  by  h+  and  k+, 
we  find,  after  a little  simplification,  that  A and  B must 
satisfy  the  two  equations 


A 


2cosa 


2B 


n sina  2 tt j n sina-nsint) 

o o o 


(28) 


A + 


2n  sinacosa 
o 


-2Bnsin0 


2 TT  j 


n sina-nsin0 
o o 


(29) 


with  solution 


A = B = 


n cosasina  nsind  -n  sin a 


”3 


nsinQ  +n  sina 
o o 


!30) 


Unfortunately , with  G(W)  given  by  (25),  it  is  implied  that 

G (W* ) ' B/(W*-W  ) and  this  introduces  an  additional  pole  into 
o r 

(21)  and  (22)  at  the  value  of  w for  which  W*  = W , i.e.  at 

o 

w =-ncos(0^  - 2a).  This  pole  requires  further  additional  poles 
in  F and  G to  satisfy  (21)  and  (22),  with  still  further  additional 
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1 

sets  of  poles  arising  in  the  same  way.  Apparently  this  series 
does  not,  as  had  at  first  been  supposed,  terminate  in  anyway, 
and  so  both  this  method  of  solution,  and  also,  apparently, 
equations (21)  and  (22),  must  be  considered  faulty.  The  reason 
for  this  is  not  known,  but  may  be  tied  up  in  the  choice  implicit 
in  (14).  However,  without  this  choice  there  seems  to  be  no 
obvious  way  to  proceed  beyond  the  setting  up  of  the  initial 
equations,  and  the  method  fails. 

2.9  CONCLUSIONS 

The  spectral  analysis  has  led  to  the  two  relations  (21) 
and  (22),  which  bear  a superficial  resemblance  to  the  more 
familiar  Wiener-Hopf  equations.  However,  an  examination  of 
their  pole  behavior  suggests  that  these  relations  are  faulty 
and  rather  reluctantly  the  method  has  had  to  be  abandoned. 


, 
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\ CRll'lQUl  Of 

,.u\nsKii's  mi  mon  of  solution  i'o  m mum  on 

riunu.i.Ms  i n\ oi v i m;  \ in  o r worn. aii  mi  i.muu:  iviikit 

MiSl'K.UM 

\ t borough  Invest  igat  ion  el'  .'a  vadsk  i i ' s method  is  made  m an  attempt  to 
obtain  solution  to  elect  romuguot ic  diffraction  problems  involving  a 
rectangular  dielectric  wedge  fl  and  f i ) infinite  metal  plate 

along  = * n / 2 . (iit  semi -infinite  metal  plate  along  4>:  n / 2 . (iiil 
perfect  magnetic  eond.uctor  along  <|>  :i / 2 and  a semi  infinite  metal  plate 
along  «+>  a In  all  the  cases  it  is  shown  that  favadskii's  method,  as 

it  is,  gives  a solution  containing  branch  cut  integrals  that  grow 
exponentially  in  the  farfield  thus  violating  the  radiation  condition, 
for  case  f i ) involving  an  infinite  metal  plate  a simple  way  of  modifying 
favadskii's  solution  is  shown  so  that  the  resulting  solution  conforms 
with  the  known  exact  solution.  Several  mod i f i cat  ions  to  favadskii's 
method  are  tried  t o obtain  the  correct  solution  to  cases  fill  and  liitl 
but  none  of  them  proved  to  be  successful . \ method  involving  a secondary 

solution,  with  branch  cut  integrals  alone,  is  shown  to  load  to  two 
coupled  integral  equations.  However  it  is  not  clear  if  a solution  to 
these  integral  equations  ousts,  finally  a unique  solution  is  given  to 
the  problem  of  illuminating  a rectangular  dielectric  wedge,  resting  on 
a semi  ml  mite  metal  plate,  such  that  there  is  no  net  diffracted  wave 
from  the  edge. 
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3.1  INTRODUCTION 

I'he  problem  of  wave  diffraction  by  a perfectly  conducting  metallic 
wedge  has  been  solved  [1,2|.  The  analogous  problem  of  diffraction  by  a 
wedge  with  impedance  faces  has  also  been  solved  [3-10].  A generalization 
of  this  solution  to  the  case  of  a dielectric  wedge,  either  free  or  resting 
on  a semi- inf inite  metal  plate,  entails  serious  mathematical  difficulties. 

\part  from  being  a classical  boundary  value  problem  the  diffraction 
of  electromagnetic  waves  by  a dielectric  wedge  is  of  particular  interest 
in  the  theory  of  dielectric  wave  guide  matching  [11,  12].  radio  propagation 
over  the  earth  [13,11],  and  in  radar,  for  the  effect  of  scattering  by 
dielectric  radomes  | 1 5 ] . An  analogous  problem  to  that  of  the  dielectric 
wedge  is  encountered  in  the  field  of  acoustics  [lo|  and  in  sei smologica 1 
situations  involving  the  behavior  of  Rayleigh  waves  at  the  boundary  between 
the  ocean  and  the  earth  [17,  18],  It  is  not  surprising,  then,  that  con- 
siderable research  effort  has  been  directed  towards  the  problem  of  dif- 
fraction by  a dielectric  wedge  [19-30], 

An  understanding  of  the  diffraction  by  a dielectric  wedge  resting 
on  a semi- inf inite  metal  plate  is  of  importance  in  assessing  the  effects 
of  dielectric  supports  for  a wave  guide  feed  illuminating  a reflector 
antenna.  A number  of  years  ago  Zavadskii  [31]  proposed  a method , which, 
he  claimed,  would  give  exact  analytic  solution  to  a class  of  two-dimensional 
wedge  diffraction  problems.  However  the  solution  obtained  by  his  method 
contains  branch  cut  integrals  which  he  does  not  evaluate.  Upon  examining 
these  branch  cut  integrals  we  found  them  to  "explode"  at  infinity  in 


complete  violation  of  t ho  radiation  condition. 

Several  attempts  were  made  to  modify  his  approach  so  as  to  remove 
this  drawback  hut  did  not  meet  with  success.  In  this  report,  after  chine 
a brief  summary  of  Zavadskii's  formulation  we  proceed  to  discuss  various 
modifications  that  were  tried  to  circumvent  the  failure  of  his  approach, 
he  also  present  the  solution  to  what  we  call,  a "qua s i - 1 r i v i a 1 " problem, 
that  ot  illuminating,  a rectangular  dielectric  wedge  resting  on  a semi- 
infinite metal  plate,  such  that  there  is  no  diffracted  wave  from  the  edge. 
Though  most  of  our  discussion  will  concentrate  on  the  problem  of  diffrae- 
t ion  by  a rectangular  dielectric  wedge  resting  on  a perfectly  conducting 
metal  plate,  we  will  also  touch  up  on  other  related  problems  to  which 
Zavadskii's  method  might  be  applicable. 

It  is  well  known  that  the  problem  ot'  diffraction  in  a wedge  with  per- 
fectly conducting  faces  was  solved  a 1 eng  time  ago  by  Sommerfeld  1-']. 
Malyuchinets  (o-7|  lias  proposed  and  developed  a method  of  solving  dif- 
fract ion  problems  in  angular  regions,  with  application  to  a wedge  of 
arbitrary  apex  angle  with  ideal  impedance  faces,  to  sectored  media  repre- 
senting a system  of  wedges  with  a common  edge  and  common  faces.  This  method 
is  based  on  the  representation  ot'  the  field  in  a dielectric  medium  by  the 
Sommerfeld  integral  and  it  reduces  the  diffraction  problem  to  functional 
equations  for  the  integrands.  But  solutions  to  these  functional  equations 
could  be  obtained  only  in  very  special  cases.  I'o  circumvent  this  dif- 
ficulty Zavadskii  introduces  the  generalized,  two-sided  laplace  transform 
which  we  denote  as  the  t-transform.  Upon  application  of  the  t -transform 
to  the  functional  equations,  and  after  some  algebraic  manipulation  one 
obtains  a functional  equation,  with  periodic  coefficients  in  the  trans- 
form domain,  which  is  amenable  to  algebraic  solution.  In  the  next  section 
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wo  give  a brief  summary  of  lavadskii's  formulation  and  in  the  subsequent 
tlttoo  sections  we  proceed  to  discuss  the  various  approaches  that  were 
tried  to  obtain  a correct  solution  to  the  diffraction  problems  involving 
t he  fo 1 1 ow i ng  g comet r i es : 

i.  Rectangular  dielectric  wedge  resting  on  an  infinite  metal  plate. 

ii.  Rectangular  dielectric  wedge  resting  on  a semi  - infinite  metal 
plate. 

iii.  \ mixed  boundary  value  problem  involving  a rectangular  dielectric 
wedge  C 0>d>>-nr/2 1 with  a metal  plate  along  = n/2  and  a 
perfect  magnetic  conductor  along  ^ = - tt/2. 

conclusions  are  given  in  section  5.(i. 
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o.d.l  Statement  of  the  problem 

We  begin  with  the  following  problem  in  cylindrical  co-ordinates. 
Consider  two  homogeneous  sectoral  media  which  have  a common  edge  IP  = 01 
and  a common  face  = 01,  each  occupying  respective  wedges  with  apex 
angles  ^ and  as  shown  in  Fig.  3.1.  The  remaining  outer  faces  of 

the  wedge  are  resting  against  perfectly  conducting  metal  plates.  l'o 


Figure  3. 1 .Geometry  of  dielectric  wedge  diffraction  problem. 


simplify  matters  we  assume  that  one  of  the  media  is  free  space  and  the 
second  medium  has  a magnetic  permeability  equal  to  that  of  free-space  and 
a refractive  index  of  n l • hater  it  will  be  necessary  to  restrict 
iJ'j  to  be  a multiple  of  ' but  at  present  we  will  assume  s',  to  be 
arbitrary.  We  assume  that  a monochromatic,  r. -polar i zed  plane  wave  of 
unit  magnitude  is  incident  on  the  dielectric  wedge  at  an  angle  ii'(1  with 
Then  the  diffracted  electric  field  will  also  be  polarized 
in  the  ^-direction  and  we  have  a two-dimensional  scalar  problem. 

We  use  the  scalar  functions  F (p , d1 1 and  F^lP,^)  where  p = kr 
to  represent  the  total  electric  fields  in  the  free  space  and  the  dielectric 
medium  respectively.  The  fields  F and  1 ; must  satisfy  the  wave 
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equations  in  fhoir  respective  regions,  the  boundary  coiniitions  .it 
A ~ l’.  4>0>  and  -A  .the  radiation  condition,  and  the  edge  condition.  V 

mathematical  description  of  these  requi rements  is  given  below. 


a)  The  wave  equation: 

U ♦ C)E(P,A)  = 0;  P<A<A„ 

(A  + k'nOl^  (P,A)  = 0;  -A.CsO 


(la) 

(lb) 


when 


3“ 

3r' 


1 3 
r T 


i y 


■V  : k = u»/< 


1 •»- 


( lc) 


is  the  angular  frequency  and  c is  the  velocitx  of  light  in  free  space. 


b)  I'he  boundar\  conditions  at  A = 0,  A , and  -A  - 

o 1 

l-  (P  »A  ) = 0 

o 

= 0 

E(P,0)  = Ej (0,0) 

E ’ (p , o)  = r! (o,o) 

where 

3 3 

1 '(PCI  = E(p,a)  ; and  l:j(P.A)  = ^ 1-  j (o  ,A ) 

1 he  conditions  (da)  and  (dbl  follow  from  the  requirement  that  the 
tangential  electric  field  be  zero  on  the  met  ll  plate  and  (del  and  (dd) 
follow  from  the  requirement  that  the  tangential  field  and  the  normal 
derivative  ot  the  field  be  continuous  across  the  boundary. 

c)  The  radiation  condition: 

■ im  r^1-  - jkf.j  = 0;  for  0 «J>«^  (3a) 

r -*«> 

- jknlij ) = 0,  for  -A(<^<0  (3b) 

1* 


him 


rh 

1 v 3r 


Ida) 

CM 

CO 

(dd) 

Co) 
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where 

F(P,4>)  = F(P,4>)  - l:0 LP ,0 ) 1 3c  1 

and  !:^(.',.>)  = exp[- ] is  the  incident  field,  and  the  time 
reference,  following  Zavadskii,  is  taken  as  e 

Physically  the  conditions  (3a)  and  (3b)  mean  that  all  the  reflected, 
refracted,  and  diffracted  fields  must  be  radially  outgoing  at  infinity. 

A detailed  discussion  of  radiation  condition  may  he  found  in  reference  [32|. 

d)  The  edge  condition  [ 1 J , 33-55]: 

The  edge  condition  requires  that  the  electrical  and  magnetic  energy 
stored  in  any  finite  neighborhood  of  the  edge  must  be  finite;  that  is, 

| { e | r: | 2 + V1 1 h | 2 } d v - o (4) 

V 

as  the  volume  v contracts  to  the  neighborhood  of  the  edge.  For  a smooth 
edge,  which  may  be  regarded  as  locally  straight,  the  differential  volume 
in  (4)  is  dv  = rdrdv'dc.  Then  condition  (4)  states  that  in  the  neighborhood 
of  the  edge,  none  of  the  field  components  of  (I'.ll)  should  grow  more 
rapidly  than  e with  t>0  as  r*0. 


3,  Z.J  Sommerfeld  integrals 

To  obtain  a solution  satisfying  (la)  - (Jd)  lavanskii  begins  with  the 
following  representation  for  E and  T ^ . 


E(P,<J>) 


:ttj 


. -jpcos(o-4>)  , 
s(a)c  da; 


• -o 


( 5a  1 


, . . 1 [ - ipncos  (<f>- . ) . . . 

P..M  = — - c - s,  ( .)  d -4'.v-4" 


-71  1 


( 5b ) 


where  \ is  the  Sommert'eld  contour  of  integration  us  shown  in 
1’he  shaded  portions  in  the  a- plane  represent  the  regions  where  the  real 
part  of  the  exponent  of  t he  integrand  in  teal  is  negative  ami  ensures 
convergence  as  a goes  to  infinity.  The  path  of  integration  'i+ls  shifts 
as  is  varied.  as  given  In  i5a)  clearly  satisfies  the  wave 

equation  since  it  represents  an  infinite  sum  of  plane  waves  each  of  which 
sutistT  the  wave  equation.  Since  the  end  points  of  the  contour  1 + 4>  lie 
in  shaded  regions  an  infinitesimal  displacement  of  the  contour  does  not 
change  the  value  of  the  integral,  which  implies  that 


£ Mp.«  - 4.  I I 4 0-j^os(»-a)1  A 
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Vhe  contour  ” is  so  chosen  that  under  the  transformat 


:,(.al  = cos 


- 1 cos  a 


with  t he  choice  of  the  branch  such  that 


, lei)  = - r.f-a) 


l,  (ocMr)  = c (a)  + it 


> transforms  to  V in  the  r.-plane. 
In  the  a.piane  I' ^ is  given  by 


fjlp.s'l  = S|  (r.  e >at  (a)  da 


SUM 


u>0 


d*  i.t ' 
cux 


> n“  cos"  i 


(SM 


Iho  !'  ranch  points  .it  * tees  'n'm  , ::i  I J,.,,  arc  joined  pair  wise 
a>  shown  in  i ig..\d,  and  the  square  root  sign  is  so  chosen  that  the  t'ollowin 
relat  ion  holds . 


T (.Ot)  = T(-Of)  (Sc) 

B>  different  iat  ■ ("cl  we  note  that 

r (cx+tt)  = t (aj  (Sd) 

and  (Sc  I is  consistent  with  (~b)  . 

Hie  relations  (~a)-(7e)  ensure  that  the  mapping  from  the  o-plane  to 
the  s-plane  is  one-to-one. 

lhe  normal  derivatives  of  1 ;i  r are  obtained  b\  differentiating 
(sal  and  (5b)  with  respect  to  0 . 


„ mm 


-.ip]  sin(o-4>)  s(ol  e'-'^'V°Ma  “Via 
1 +0 


dtt  i 


- -.ip  /si no  s (o+pi c 


- -jpfnsinu-vMsl(rje-j‘'ncos^\lr. 


r +4" 


. „ r . , - ipneos'.  , 

-iPJnsin:.s1(..+pH'  ' d;. 
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he  have  here  made  use  ol'  the  relation  t(a)  = sina/nsinC. 

Using  the  expressions  (5a),  (5b),  (9)  and  (10)  for  E,  E, , and 

1 <3 

3Ei  ... 

respectively  we  obtain  the  following  set  of  functional 

equations  for  s and  s,  so  as  to  satisfy  (2a)-(2d). 


s(a+<p0)  = s(-a+4>0) 


s1(C(a)-lp1)  = s^-Ua)-^) 


s(a)-s(-a)  = t (oc)  (s1[C(a)]-s1[-5(ot)]} 
s(a)+s(-a)  = s1[?(a)]+s1[-C(a)] 


(11a) 

(lib) 

(He) 

(lid) 


3.2.5  Generalized  Two-sided  Laplace  Transform 

In  order  to  solve  (I  la- 1 Id)  -avadskii  inti'oduces  the  following  trans- 
form relationship. 


0(a)  = 


t (a,  p)  e 


"UPdB  ^o(a)  = t(a,3)j  (12) 


which  represents  a function  of  one  argument  in  terms  of  a function  of 
two  arguments.  If  t(a,0)  is  not  a function  of  3 then  the  above  trans- 
form goes  over  to  ordinary  two-sided  Laplace  transform  [36].  It  should 
be  noted  that  given  o(a)  the  function  t(a,0)  is  not  unique  since  to 
any  given  t(a,0)  satisfying  (12)  one  could  add  any  other  function 


t^(a,3)  with  the  property, 


to(ct,3)e'a0dg 


(13) 


In  the  case  when  the  function  a(a)  has  a pole,  at  the  point  a=i^,  with 
unit  principle  part,  in  the  strip  Re<KRevo<Re(a *J>y) , we  then  have 


oo 


o(a+4>0) 


0^()  *ii  3-<t>03 

| t [a  » p ) + c ]o  <J3  ( 1 4 j 

-oo 


The  second  term  In  the  brackets  corresponds  to  the  pole  of  the 
function  at  the  point  a = ij;  . The  convenience  of  the  representation  (Id) 
lies  in  the  fact  that  in  certain  problems  the  function  t(«,l3)  may  be 
regarded,  without  contradicting  the  conditions  imposed  on  the  problem, 
as  periodic  in  the  argument  a,  whereas  the  function  0(a)  is  not  periodic. 
An  example  is  the  case  when  0(a)  satisfies  the  functional  equation 


o(a)  = q(a)o  (a+(j>Q)  + f(a) 


(15a) 


where  q(a)  = q(a+rr)  and  f(a)  is  representable  in  the  form 


f (a)  = F(0)e 

-oo 

-a0  ,0  . 

dp,  <p()  = constant 

(15b) 

T 

ak  i ng 

t he  t runs  form  of  (1 5 

a)  we  obtain 

A)3 

t (a , 3 ) = q (a)  e t(a+<f>0,6)  + P(8) 

(l(h 

whose 

solution  may  be  written  in 

the  form 

t(a,8)  = F 

-<J)  8 

C3) L [q(a)e  0 ] 

(17) 

where 

I-(b  (p(a) ) " 1 
f() 

°°  k 

+ l l p Ccx+v^h  ) 
k=0  v=() 

(18a) 

p(a)  = q(a) 

C 

(18b) 

S i nc  e 

q (a) 

= q(a+iT)  , it  follows 

that  | p (a) 1 is  also 

periodic  with 

period 

IT  . 

In  special  cases,  assuming  that  the  series 

(18a)  converges. 

bl-’ 


wo  have, 


I [p(a) | I/ll *p(a) | 


V y"  mv 

I + 11  Plot*  ■') 

'll  V l' 

5 i 

1 p (0t  + IUX  IT  1 

V 0 


m = 1 , 


« 1 


Now  wo  turn  to  t ho  t'tinot  i on;i  1 equal  ions  (llal  - (lid)  and  invoko  tho 
following  roprosontat  ions  for  s(a)  and  s (ix)  : s(cx)  i t tot . k> ) . 
s | (•'.  (>.x) ) ' tj(,ix,3),  witli  tho  ros  t r i o t i on  that  if  is  an  integral  multiple 
of  n/2.  Continuing  tho  roprosontat  i on  s(ol  = t(cx,|3)  beyond  tho  polo 


at  ot  1 , wo  obtain 


Sloi-s^,'  1 1 ux+4>0,tn  + o I o do 


Since  t (>xl  r (a+n ) , tho  functions  t(a,S),  t (ix.fl)  have  a period  n 

with  respect  to  a.  With  this  fact  in  mind  we  obtain  the  following  system 

ot  inhomogeneous  functional  equations  for  t(o,b)  and  t (ix,b): 

bij't,  -O  B>j'0  ^ b 

l 1 l'xuji(),d)  i o | e |t(-<xn|>(1.  dl  f c | e C --a  1 


t t (ix.n  ) t l -,x , die 


t bx.b)  tl-o.-B)  t (a)  1 1 . (iX.d)  - tjl-a.-dl 


l -’-'hi 


t-'.'cl 


i (o.d)  *■  1 1 -u,  b ) 


l | l'X.l'1  * t j l-o.-b) 


which  can  bo  reduced  to 


--V()o  -<!>  b 

t(.o,p)  [t (.ot+jy.b )e  + Jo  shbu  | 


l-i  la)  tlujy? 
1+t (a)th<{>  3 


1 ho  solution  to  C — o ) may  bo  written  in  the  form 

V 

t[oi,bl  Jo  slUiJ<  -v|>  ) 

ho  may  now  obtain  s (pV)  by  inteyratiny  (Jl)  . 


H 1-  t (cdtlnfy! 

1 

x 1 

[ 'J""° 

i~f(aTt  lu,  t b J 

-» 

J 

IJD 


s (.a)  = 


t (ot , b ) o 1 do 


125) 


lo  find  the  function  s y.  (et) ) , which  in  turn  defines  s ^ ti-n ) implicit  lv, 
we  may  either  find  tyx.bl  through  iJJc).  (JJd)  and  (Jl)  and  intoyran 
or  substitute  for  stall  in  tllcl  and  tlld)  and  eliminate  s (yul ) . 

In  principle  one  could  invert  (Jll  for  any  arbitrary  with  41 

beiny  a multiple  of  n/2.  However  when  y is  not  a multiple  of  n/J,  i 
ia>  not  be  possible  t <■>  obtain  a closed  form  expression  for  s(oO  throuyh 
125). 

J.l  A note  on  the  solution  to  functional  equations 

1 von  thouyh  stal  and  sy-to)),  as  obtained  by  the  above  procedure, 
will  satisfy  (lla)-llld)  there  may  be  other  solutions  to  this  system 
ot  equations.  There  are  two  reasons  for  this.  The  first  is  that  t (yt.b) 
as  yiven  by  (Jit  may  not  be  the  only  solution  to  tJJa)-(JJd).  The  second 
and  the  more  important  reason  is  that  as  we  mentioned  earl ier, yiven  sp\\ 
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the  transform  t(,ix,|i')  is  not  unique.  lienee  the  oqns  l J Ja ) - (JJd)  do 
not  represent  a uni  quo  transform  of  t ho  functional  equations  1 1 1 a ) - ( 1 1 d) . 
llonco  the  solution  to  s(ot)  anil  Sj('.(a))  as  ohtainoil  above  Joes  not 
constitute  a unique  solution  to  the  problem  but  it  is  only  a particular 
solution  to  the  functional  equations  (11a)  (lid).  To  check  the  correctness 
and  the  validity  of  the  solution  one  must  evaluate  1 (,  .s')  and  1 ' (,  .^l 
and  verify  if  the  fields  so  obtained  satisfy  the  radiation  condition  and 
the  edge  condition.  Hut  as  we  shall  show,  through  specific  examples  later, 
the  fields  obtained  through  ..a vadsk i i ’ s method  diverge  at  infinity  and 
hence  fail  to  meet  the  radiation  condition.  When  i|>  is  a multiple  of 
tt/J  it  is  possible  to  obtain  the  most  general  solut ion  to  (lla)-(lldl  with 
out  taking  recourse  to  t-transform.  However  any  solution  to  this  system 
of  equations  is  found  to  result  in  either  an  undesired  incoming  plane 
wave  or  a branch  cut  integral  with  exponential  growth  at  infinity.  both 
or  either  of  these  conditions  constitute  a violation  of  the  radiation 
condition.  Since  the  system  of  oqns  (lla)-(lld)  and  the  representation 
of  the  fields  as  a Sommer  fold  integral  over  the  path  y are  the  ke\ 
steps  in  favadskii's  method  we  concluded  that  his  method  tails  to  give 
the  correct  solution  to  the  diffraction  problem.  We  demonstrate  this 
through  a few  specific  examples. 


3.3  RHCTANTUl.AU  HIl.I.LCTRU'  IVlilK’.r  RHST1  NC  ON 
AN  I NF I N IT!:  MITAL  PLATL 

3.3.1  Zavadskii' s solution 

Lot  us  consider  the  case  when  <’  ^ ^ = tr/d,  for  which  the  exact 

solution  is  well  known  through  geometric  o]itics. 

In  this  case  (d-1)  simplifies  to 

2sh3(i^-  t)  ‘ ch^o-r(ct) sh  -,-3} 

{ l+i(al } s h it 3 (del 

which  upon  integrating  results  in  the  following  expression  for  s(o) 


At  this  point  Zavadskii  stops  with  a somewhat  bland  statement  which 
we  quote  here.  "The  poles  of  the  function  s(a)  situated  in  the  strip 
-tr<a<TT  for  a = ' i Ji  {witli  principle  parts  equal  to  1 and  1 1 -T (iji  1 ]/ 
[l+x(i|i  )),  respectively}  and  for  tx  = (it-tji.j)  {with  principle  parts  -1, 
and  - | 1 -r  (i|j() ) ] / | 1 + 1 (i/;  } } make  it  possible  to  compute  the  Sommerfold 
integral  (5a)  and  to  obtain  the  field  in  the  wedge  (0,ir/d)  in  the  form  of 
a sum  of  four  plane  waves.  In  this  case  there  is  no  cylindrical  wave 
radiated  by  t he  edge  of  the  wedge." 

However  Zavadskii  overlooks  a very  important  factor.  lo  evaluate  the 
integral  (5a)  one  has  to  close  the  contour  y in  order  to  pick  up  the 
poles  lying  on  the  real  axis.  One  would  normally  do  this  by  adding  two 
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paths  of  steepest  descent  D ,D  which  are  deformed  around  the  branch 
cuts,  due  to  x (,oO  , as  shown  in  Fig.3.3.  In  this  particular  case  D ,1) 
are  spaced  apart  by  exactly  Jtt  and  s (a ) being  periodic  with  period  Jit 
the  integrals  on  D and  P cancel  each  other.  Then  we  are  left  with 
the  branch  cut  integrals  (along  the  paths  B , B^  as  shown  in  big.  1) 
in  addition  to  the  residues  from  the  poles  at  a = ±ip  and  ct= 

As  we  see  from  Fig.3.1  the  upper  portion  of  B f and  t he  lower  portion  of 
B lie  in  the  unshaded  region  where  the  exponent  of  the  Sommerfeld 
integrand  lias  a positive  real  part.  Tluis  both  of  these  branch  cut  integrals 
grow  exponentially  as  p-*00  and  fail  to  meet  the  radiation  condition.  In 
Appendix  l we  show  that  the  branch  cut  integrals  over  the  paths  Bj  and 
B,  do  not  cancel  each  other.  Thus  we  conclude  that  Zavadskii's  solution 
is  non-physical  and  is  of  no  practical  use.**  We  know  that  for  this 
particular  geometry  the  exact  solution  is  just  the  sum  of  four  plane  waves 
as  obtained  by  the  geometric  optics  method.  As  we  noted  earlier  the  solution 
obtained  through  the  t-transform  is  one  particular  solution  to  the  system 
of  equations  (lla)-(lld),  but  not  necessarily  the  correct  one.  Any  correct 
solution  to  our  problem  in  addition  to  satisfying  these  equations  must 
also  satisfy  the  radiation  and  the  edge  condition.  It  turns  out  that  there 
is  in  fact  a solution  to  (lla)-(lld)  which  meets  these  requirements.  We 
proceed  to  obtain  such  a solution  as  follows. 

3.3.2  Correct  solution 

Let  s(al  and  s ^ (a) ) meet  the  following  requirements. 


* * 


This  was  pointed  out  to  Zavadskii  in  a personal  communication  to  which 
there  was  no  response. 


Contours  of  integration  > , > , n , ami  I'  in  the  plain 
o1‘  t ho  complex  variable  a. 


;>-i9 


il  s(a  + ?r)  = s(-a) 
ii)  sia)  = s(a  + 2?r) 

iii}  V5(«j)  - 1 sl  “ 


S (. 01  ) + S [ ) 


Then  si'll  and  s iC.(eil)  will  also  satisfx  (llal-(lld)  with 
0()  - = i/J.  \ny  odd  function  with  period  2^  or  any  even  function 

with  period  it  will  satisfy  (2Sa)  and  (2Sb).  Thus  one  lias  an  infinite 
number  of  solutions  to  ( 1 la  1 - ( 1 Id) . If  we  impose  the  condition  that  the 
result  ini’  field  must  have  only  one  incoming  plane  wave  at  <{•  we 

narrow  down  the  solutions  to  two.  One  solution  is  that  given  In  (2~). 

But,  as  already  noted,  it  fails  to  meet  the  radiation  condition,  and  hence 
is  not  the  desired  solution.  The  second  solution  which  unique!)  meets 
all  the  requirements  is  simply  obtained  In  replacing  T(o)  with  T C ^ 1 in 
(27) . Titus 

cosV(,  1 -T  cos  Hi 

S(.u)  = — : ! T r : T- (29a) 

sum-  s i n,. ^ l + ti^^l  sina+sunji^  v 


Substituting  (29a)  in  (28c)  we  obtain  the  following  expression  for 
SjUln) ) • 

2rU'0lcosy0  sin^/rU^  + sina/x(a) 

»|(5(o))  = i+t  (ij  j 71  l 

i'  sin  a-  sm 

o ; -o'  * nsn';-  (29b) 

1 M^V  cos“t  - n“cosT, 

0 


It  immediately  follows  that  s.pt)  is  given  In 


(•*)  = 


2T(*0)cosV0  sinvt)/T(v  ) ♦ ns  inn 

- ■ - C 

1 r I^'q)  cos*>  - n“cos“« 


2x(»0)cos60  sing  ♦ sin90  2Tty0)cos  qq  t 

* cos',1  -cos',-,  1 * sm  i -.in  ' 
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whore 


cose  = cos-,  ,/n 
0 0 


l-h'd) 


It  is  now  cas>  to  verif\  that  s(0}  and  s (g  (.x) 1 , as  given  by  (2da,bl 
do  indeed  satisfy  (llal-(lld).  Since  the  onl\  singularities  of  s(al 
•lie  the  poles  at  o = .*  i;^  and  ot  = we  have  no  branch  cut  integrals 

to  ev a lu.it e and  the  electric  field  1 (P.P!  is  simple  given  bv  the  follow- 
ing sum  four  plane  waves. 


-.iPcos  (v'-vl,  1 jpeos (p+y ,,1 

1 l e - e 


■jpeos  (P  + , ^1 


jpcosW-^l-j 


(301 


Mso  since  s.(o0  does  not  have  any  branch  cut  singularities,  the  electric 
tield  inside  the  dielectric,  1’.j(Pj40  is  given  by  the  sum  of  two  plane 
waves,  due  to  poles  at  a = 0 . ir  -0 

\fter  evaluating  the  residues  we  find  the  following  expression 
for  I;  (o ,*): 


-t(P(1l  , -jpncos(<}>-0  ) jpncos(P+0  ll 

'h1^1  "fnv  Y I tM) 

Hie  tield  expressions  1:  ;i  p as  given  by  (301  and  (311  correspond 
exactly  to  the  geometric  optical  solution  to  our  problem.  Thus  we  have 
shown  a way  to  overcome  the  difficult'’  in  finding  the  solution  to  this 
problem  through  -avadskii's  method.  Unfortunately , however,  the  same 
kind  ot  approach  does  not  work  to  other  non-trivial  geometries  as  we  will 


demon st rat  e. 


RtA'l  Vxi'.Ul.AR  I'll' Id  i'TUH'  IMlV.f  RISl'lNC  0\  \ 
3!  Ml- l\!  l\l l!  Ml  ! \1  l'l.AIT 


1.1  lav.ulskii's  solution: 


Consider  the  ease  where  a semi -infinite  metallic  plate 


is  resting 


against  a 

rectangular  di 

electric  wevlge  as 

correspoiu 

.V  i 

II 

•/. 

and  = n/2  in 

Hence  all 

the  equations 

t h rough 

t.  2 J A hoUl 

for  this  probl 

em.  Sub- 

st itut  mg 

the  values  of 

s''  . and 

0 

C j we  obtain,  after  some  simplif- 
ication, the  following  expression 
for  the  function  t|..\,31* 


shown  in  l;ig,\5. 


l'his  situation 


Fig.3.1 . 


figure  3.  5.  Rectangular  dielectric  wedgi 
resting  on  a perfectly  eon- 
ducting  semi  - inf ini te  plate 


tpx.el  = 2sh(V  - V:  (ch  It'  - T(o0  sh  -JoJ 


( 1+T (a) "J -.h;  * [l-i  [col shit 3 (321 


As  shown  m Appendix  11.  eqn.  (321  leads  to  r lie  following  expression 


for  the  functions  slul  and  s^  i.:.  l.\)  1 : 


s(a)  \ Uv./o-nl  + \ lu  a-:-0 
l 1 „ O 


* A I U 0-o- il  + At  -) 


l 33a  1 


( 3 3b ' 


where  A and 


are  vie  fined  bv 


\ 


cos  t 


(..'.V  ' 


mu  1 0~  I 


l'' 


! sh  w d\ 

; sh  v iehfiy+;\) 


I , . - \ (.  r - 1 ) 
. .sin 


- Cl  Ms  in  . 


I 1 - \ “ 1 s i n \ 


(Pad 


Vhc  nolcs  of  Seal  are  located  at  a ' - .-in:!  where  m is  an  integer 
■ 0 


lo  evaluate  the  Sommer t*e Id  integral  iSal  we  close  the  contour  \ In 
means  of  two  steepest  descent  paths  P and  D as  shown  iit  fig.  a. 
Then  l'l.p.v’)  is  given  by  the  sum  of  the  following  three  terms. 


hld.O)  1 (0.0)  * + fjCO.M 


l-'ll 


where 


yo.c) 


- ^ Residues 


1 00  a 


l , ipcos  ia-uM  , 

, i soil  e ■ da  lKa<n 

dn  i J 


yp-*> 


B +R 
0 n 


i 2TT  i 


. -jOcosia-01 
s(.a)e  da 


a!T 


l a oh 


>'•!  fR  >!! 


h/0.0)  = 


-1_ 
-hi  I 


s la)  e 


- .ipcos  ia-^) 


d a 


l arc 


n +n 


J.l.l  pole  contribution 

The  steepest  descent  paths  P 


and  P intersect  the  real  axis  of 
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a at  - 'I  + and  ' + <j>  respect  ivoly . Accordingly  the  combination  of 
poles  enclosed  by  the  contour  changes  as  is  varied.  After  some  algebraic 
manipulation  (.Appendix  111)  s (ex)  may  he  put  into  the  following  form 
which  is  more  convenient  for  the  evaluation  of  the  pole  solutii 
t [2A”+A-  l)  (.2s inocosC'^l  - (l+A-2.V)sin2’.( 


:ion  L (p,d>): 
l1 


cos2x\  - cos2v 


(l -A  1 (cos 2a-  co$2v^)  l 


1 . , /A'1  5t\  j'xt  5t\ 

<sm.v0cos^—  - — jcos^~  - 3-) 

1 

. (AA  3t\  . /ext  5t\; 

+ - j-Jsin^—  - rj. 

-A 


- f)“s  (g 


(3b) 


The  residue  R(o^)  at  any  pole  cc  = is  then  given  by 

2s  ina  cosV  ,1 


l 


A Sin-Vt1 

/ > \ 

(1+V2Aoj  - 

; 1 n_o  . 


where 


* Jsln-vo  (Vo  5t0\  (% 

, 0,3  , ) 


,a"ro  -'A 


a -a-  i 


A)  ' 


and  A0  = A(i^) 


The  residues  of  poles,  that  are  of  interest  to  us  are  given  by 


(5“a) 

(57b) 
( 3 7c ) 


i 


R«V  =l 

Rt37r-%^  = -1  (38) 


1 ho  location  of  poles,  along  with  the  residues,  are  shown  in  Fig.  6 for 
three  different  ranges  of  the  incident  angle.  The  points  of  intersection 
of  D_  and  D+  with  the  real  axis  of  oi-plane  are  also  shown.  Bv 
looking  at  Fig. 36  we  can  immediately  write  down  the  pole  solution  E(p,4>) 
for  different  ranges  of  cf>  and  as  given  below. 


a)  0 L \fr( 


4 A o - I 


L.  JL 
2 


X 


-H  *0k- 

2A0  ! i 


I 


k. 

A 


-1 


a .!') 


i i 0*  t^'  , 


I ..IP  .1’  I 


\ 


''  *0 

> I. 


4', 


, i I tp0-  , .1 


III) 


M"S)V 


.111 


I (p.,|0  |0 


' 0 

r r 
% " '(> 


will'  IV 


I’ 


ii’i'os  (.;>  < 1 ) 


; , , hJ 

i)  ii 

3 IT 


, 'i  'A»p  w, 

1 

iM»<n  •' 

1 

/ p 

n -tfi  * 

\ V0 

1 0 

p -p 

2 n - tf; 

( % n % 

O 

, , .1 II 

• :l  V 2 


( aPa  1 


l aPb  ) 


faPe ) 


( ,'iM  1 


Win'll  llu'  ineidont  wave  is  in  tin'  first  i|uadrant  wo  observe  tho  i no  i ilont 

wn  vo  .uni  tho  roflootoil  wavo  from  tho  ilioloctric  boundary,  tint  il  i(i 

.i|'|M'o.ii  lios  i In-  valuo  n alter  wliioli  tho  rot  loo  toil  wavo  il  i sappoa  rs  . 

Klion  i)'  o woods  tho  valuo  ii  i iji  ( tho  i no  i ilont  wavo  also  disappears  ami 

wo  a ro  in  tho  dark  rop ion.  W'hon  ip  is  in  tho  soooml  quadrant  tho  i no  i ilont 

wavo  exists  in  tho  ontiro  ranpe  0-  .)-•  .mi/ .1  and  illuminates  both  tho  faros 
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of  the  wedge , thus  Riving  r iso  to  two  reflected  waves.  The  reflected  wave 


t rom  the  dielectric  boundary  exists  in  the  range  0<4><tt-vJ>  } and  the 
reflected  wave  from  the  metal  plate  exists  in  the  range  2n-iji  <$<3Tr/ 2 . 

''hen  ;|‘t)  is  in  the  third  quadrant  we  are  in  the  dark  region  until 
approaches  the  value  i^-tt  after  which  we  just  see  the  incident  wave, 
linen  T exceeds  the  value  the  reflected  wave  from  the  metal 

plate  also  comes  into  the  range  of  ihservation.  lints  the  pole  solution 
is  consistent  with  the  geometrical  optics  theory. 
a.I.l.J  Piffracted  field 

The  integrals  over  the  contours  P and  l>+  represent  the  diffracted 
waves  due  to  the  edge.  We  choose  the  paths  l)  and  P such  that  the 
following  conditions  are  met  for  ct  = p+  iv  ; P_  ,P  . 


Re{cos  } = cos  (.p -if ) chv  = -1 

Imfcos } -sin(b-i)>)  shv  m) 


l hen  u and  v are  defined  through  tlie  relations 


p = -gd  (,V)  + >J)  *-  ir 
and  sinfp-vf'j  = + t h v 


where  gd(,.x)  is  the  Hudormann  function  given  by 

gd(x)  = cos  (l/chx)  (He) 

We  should  note  that  P and  P as  defined  through  iUa.li)  mav 
hut  the  branch  cuts,  for  certain  ranges  of  n*>,  in  which  case  we  deform 
P and  P around  these  branch  cuts  as  shown  in  Tigo.a,  and  the  result 
ing  branch  cut  and  end  point  integrals  must  he  evaluated  separately. 
However  for  large  values  of  p the  contribution  from  these  integrals  i> 
negligible  compared  to  the  dominant  contribution  from  the  vicmin  of 
the  saddle  points  at  a = (p  * it . Now  wo  proceed  to  evaluate  these  saddli 
point  contributions. 
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Let;  Lj_(p,(|))  and  Ed*  Lp , 4> ) denote  the  integrals  over  0 and  H 
respectively.  Then 

1 ( . . > -jpCOS  (yt-dj) 

stale  1 da 


LJ+lp.v^  = 2rrJ  J 


J 

-n  i 


s ta)  c”  jp  ^cos  (u-4>)ehv- j s i n (p-4>)  shv  | / dp  + A 

\dv  " J;  v 

~co 

= !A-|AfoH'jP'PthUahY-l/cl,v*jlJv 


(42a) 


Now  we  approximate  s(a)  (-1/chv+j)  by  its  value  at  the  saddle  point 

•) 

a =tt  +()) + j 0 and  thvshv  by  v“  to  obtain 


'•'d+Cp,‘M  = 3i7j  s(ir^)(-l+j)  [ eJP'°V  dv 

u 

= s(TT+<f))eJ  (P+  4 //3rfp 


(12b) 


Similarly  we  obtain  H as 

d - 


l.j_(P,<J>)  = - s(-ir+s‘))e-' (-p+  dV/dTtp 


and  taking  negative  of  the  sum  of  (42b)  and  (42c)  we  obtai 


f , (n  ,4>  1 - - { s (it  +<)> ) - s ( -it+4>  ) } e*’  (p+  41/>/“"P 


(42c) 


(1 

- s 111241s  in' 


“S1Utl  . . „ j 

— — 7 sin 2ty  cosl 

A )(cos2<f-cos2tJi  ) I ‘ 

1 o ( 


Vi  5,i\  /+>  1 


it  " 2 r"'l  it 


Vi  3ti\  1**1  3t  A 

IT  2 /COS\  it  2 J 


(1-Aj  ) 


(Vi 

ini 

3tl) 

1 1 

2 J 

s 1 11  1 

\ n 

- — I 

J(P*  t) 


■ > 2np 


(43a) 


where 


t,  = t(*) 
and  Aj  = A (4>) 


(43b) 
1 4 3c) 
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If  the  refractive  Index  of  the  dielectric  medium  is  made  to  approach 
unity,  thus  reducing  it  to  free  space,  then  (43a)  simplifies  to 


j(p+  j)  ; 


l 


ijipail  ■ k w.*  • - < «,n. 

n=l  2/2tto  I i)  0 

:sin  — — — cos  — r — 


(44) 


which,  with  proper  interpretation  of  <J)  and  corresponds  to  Sonunerfeld' s 

solution  [J.  | to  plane  wave  diffraction  by  a perfectly  conducting  semi- 
infinite plate. 

In  passing  we  should  note  that  1:,(P,4')  as  given  In-  (43a)  is  onlv 

U 

good  when  there  are  no  poles,  of  s(a),  in  the  vicinity  of  the  saddle 
points  a = <j>+ir  . 

3.4. 1.3  Branch  cut  integrals 

The  contribution  due  to  the  branch  cuts  enclosed  within  the  contour 

Y+D  +D  is  given  by  (35b) . If  we  try  to  evaluate  these  branch  cut 

integrals  we  find  that  the  top  halves  of  B,,  and  B,  and  the  bottom 

0 -IT 

half  of  B ^ give  rise  to  terms  which  increase  exponentially  as  p in- 
creases and  thus  these  terms  fail  to  meet  the  radiation  condition.  On 
these  branch  cut  contours  the  exponential  term  in  the  Sommer fold  integral 
(5a)  lias  a positive  real  part  which  is  proport iona l to  psh(|v|)  where  v = I mo 
It  can  be  easily  verified  that  these  exploding  branch  cut  integrals  do 
not  anihilate  each  other  and  the  presence  of  these  terms  makes  the  solu- 
tion, obtained  through  Zavadskii's  method,  physically  meaningless.  In 
the  next  section  we  describe  several  approaches  to  try  to  correct  this 


problem. 


3.4.2  Attempts  to  correct  lavadskii's  solution 
3. 4.2.1  Is  a branch  cut  free  solution  possible? 

If  we  replace  T(a)  by  T(if())  in  (.30)  the  resulting  function  whicl 
we  call  s (ct)  still  satisfies  the  equation  (11a)  with  if  ^ = 5tr/2. 


3 it,  5tt 

s(a+  — ) = §(-a+  -s-) 


thus  meeting  the  boundary  condition  at  if  = — . However  the  modified 
function  s.(C(ot))  as  given  by 


1 | „ _ i 

A l ■'■.(«))  = Trrrr  S s(a)-s(-oi)  + T(a)|s(a)  + s(-a)]>  (,4b) 

t 2T(.a)  | 


does  not  satisfy  the  boundarv  condition  at  if  - -tt/2.  That  is 


s t (?,(a) ) / (C(-a-tr) ) 


We  note  that  s(a)  as  obtained  above  is  non-periodic.  In  Appendix 
IV  we  show  that  any  branch  cut  free  solution  s(a)  to  (11a)  with 
4>()  = 3 tt/2  must  be  periodic  in  a,  with  period  2tt,  so  that  the  result- 

ing function  s^(C(a))  will  satisfy  (lib)  with  if  = tt/2.  However,  as 
we  shall  show  later,  any  such  non-trivial  periodic  solution  would  give  rise 
to  more  than  one  incoming  plane  wave  thus  violating  the  radiation  condition. 
Thus  it  seems  impossible  to  obtain  a solution  to  the  functional  equations 
(lla-lld),  with  if  = 3tt/2  and  if^  = -tt/2,  such  that  s(a)  is  free 
of  branch  points.  I'uther,  we  know  that,  such  a solution  would  not  have 
anv  diffracted  waves  since  I'  and  H are  then  spaced  hv  2ti  and  are 
in  opposite  directions.  We  know  that  for  the  problem  under  consideration 
there  must  be  diffracted  waves  from  the  edge  and  the  associated  lateral 


waves.  Normally  one  would  expect  to  obtain  lateral  waves  from  the  branch 


cut  integrals;  hence  seeking  a solution  which  is  complete!)  free  of  branch 
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cuts  seems  to  be  ;»  step  in  the  wrong  direction.  For  the  trivial  problem 
considered  in  the  earlier  section  we  knew  that  there  would  he  no  diffracted 
or  lateral  waves  and  our  attempts  to  obtain  a branch  cut  free  solution 
then  met  with  success. 

• l*-.-  is  there  a solution  with  zero  contribution  from  the  unshaded  part 
of  the  branch  cuts? 

lo  answer  this  question  one  must  look  at  the  most  general  solution 
to  (.1  la-1  Id  1 . he  have  obtained  such  a general  solution,  in  Appendix  V, 
which  is  given  by 


s(a)  = J:7r(00 


P la)  cost 


at 


en 


% 


, at  5t , 

l)nsln(T'T) 


(,48a) 


where  P^  and  P^  are  any  even  and  odd  functions  of  a with  period 
it  and  J (a)  is  any  function  of  a with  period  Jtt  and  satisfying 
the  condition 


3 JT  3TT, 

Ca+  — ) - —) 


1 4 8b) 


1 he  function  t (ct)  is  del  ined  as  in  (33c).  It  is  ease  to  see  that 
s(a)  as  given  by  (3o)  is  indeed  in  the  form  of  (18a). 

'he  periodic  part  d,((a)  in  (48a)  may  be  chosen  such  that  it  is 
tree  of  any  branch  points.  However  the  non-periodic  part  contains  branch 
points  due  to  t.  In  Appendix  VI  it  is  shown  that  for  the  branch  cut  con- 
tributions from  the  unshaded  branches  to  vanish  the  functions  P and 

ett 

Pp  must  identically  vanish.  As  we  shall  see  a little  later  any  solution 
of  the  form  dlrt(a)  gives  rise  to  more  than  one  incoming  plane  wave  which 
is  again  a violation  of  the  radiation  condition  which  requires  that  there 
be  no  incoming  waves  other  than  the  incident  wave.  thus  there  is  no 
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solution  spi)  which,  after  integration  over  the  contour  ),  would  give 
fields  satisfying  the  radiation  condition,  and  we  conclude  that  > is 
the  wrong  contour  to  start  with.  This  immediately  gives  rise  to  the  follow- 
ing questions:  Kluit  other  possible  contours  are  there  to  s t art  with, 
so  that  we  could  write  down  t he  boundary  conditions  as  a set  of  functional 
equations  for  s(ot)?  Is  it  possible  to  combine  the  solutions  obtained 
through  ) and  some  other  contours  in  such  a manner  that  the  radiation 
condition  is  uniquely  met?  We  will  make  a thorough  examination  of  these 
possibi 1 i t ies  in  later  sections,  but  before  we  do  that  we  would  like  to 
make  the  following  comments  about  the  periodic  solution  d , (m) . 

If  J,  lot)  is  free  of  singularities  then  the  resultant  fields  would 
be  identically  zero.  Any  branch  cut  singularities  involving  t (a)  would 
give  rise  to  exploding  waves  at  infinity.  Thus,  ,I,„(o)  can  only  have 
poles.  We  specifically  require  a pole  at  o = ■P  ^ corresponding  to  the 
incident  wave.  However,  because  of  the  periodicity  and  the  requirement 
1-lSb) , d ,n (a)  cannot  have  an  isolated  pole.  Any  pole  with  a residue 
'a'  at  a =ip  ^ gives  rise  to  a chain  of  poles  with  residues  as  indicated 
in  fig. 5.7a.  It'  we  introduce  a pole  with  a residue  Mi'  at  oi  -t  tins 
gives  rise  to  a different  chain  of  poles  as  shown  in  l:ig.5.?h.  Thus  the 
most  general  pole  structure  of  J,  pi)  is  as  shown  in  rig..\?c,  Since  the 
polos  at  and  -v  correspond  to  the  incident  wave  and  the  reflected 

wave  from  the  dielectric  we  require  that 


a = l 


and  b = - 1 1-t  U>())  l/[l*Ul|>0)  1 


Then  ,I,^(oi) 


is  given  by 


I , 

7 IT 


(o) 


cos  v 


0 


s in.t-  s i n.  • % 
v 0 


1 r(V  cos»Q 

1*T(*0)  ';in  sin'i'0 


1 10a  1 


Cl01>) 


l 50) 


A 
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live  resulting  solution  for  the  fields  is  given  by 


E(p,<ii) 


- j P co  s 1 P - V ' j 1 


1-TUi'  ) -jPCOS  (<£+>,  J ) 

l+r"lv  L' 


jpcos  ()) 


+ i+uv0)  0 


(.51 


But  1*  ( p , y 1 has  two  undesired  incoming  plane  waves  at  angles  tt  - e , 
and  a + y^,  and  is  not  t he  solution  to  our  problem.  However  in  the  process 
of  the  above  discussion,  we  have  found  a unique  solution  to  a 'quasi- 
trivia 1'  problem  which  wo  give  next. 


5. 4.2.5  Solution  to  a 'quasi-trivial ' problem 

Can  we  illuminate  a rectangular  dielectric  wedge,  resting  on  a semi- 
infinite metal  plate,  in  such  a manner  that  the  resulting  field  will  have 
onlv  plane  waves  and  no  diffracted  waves?  live  answer  is  ves,  and  the 
solution  is  the  following  combination  of  incident  waves. 


H. 

me 


(P,4>3 


1-T^'p) 

1+1 1'"  o1 


r 


(.521 


where  1'^  = exp  [ - jpcos  (4'-0)  j represents  an  incident  wave  of  unit  ampli- 
tude, at  an  angle  0.  The  three  plane  waves  1' 

*0 

(R  = - [ 1 -t  ly’jl  1/ [1  ♦tl'-yl  |)  and  tlve  corresponding  reflected  and  refracted 
waves  produce  continuous  fields  in  the  entire  region  iK^<  2:r  as  shown 
in  l'ig.  S,  where  the  incident  waves  are  shown  with  solid  lines  and  the 
rot  looted  and  refracted  waves  are  shown  with  dotted  lines.  The  discontin- 
uity of  lj  ) at  ip  = 7T«-v(i  is  cancelled  by  the  reflected  wave 


and  -Rl’ 


Tt«V, 


1 1 him  in.it 
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.1  roc  l 

.until  1 ,i  r .1 1 o 1 oc  t r i c 
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.uni 

rofractod  pl.ino 
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R,01=1\,  ) and  r !\o  discontinuity  of  l,(=-RP,  ) is  nullified  bv  the 

v0  J 

ref  looted  wave  K , , (_  = — R 1 ' 1 which  exists  onlv  in  t He  ranee  0 < < ■, 1 

-1  u+'*o  ' 0 

Ihe  incident  wave  1,  gives  rise  to  a reflected  wave  R_  which  nreciselv 
fills  up  the  shadow  region  of  R^  . The  sharou  region  of  T,  is  filled 
up  b.v  fjp  which  is  the  reflected  wave  corresponding  to  f . 

It  we  let  f 1,0  to  be  the  diffracted  field  due  to  an  incident  wave 
1'  , then  1 1,0  must  satisfy  the  following  functional  equation. 

f 1 -0  - f 1 rr  - ip  1 - Rflft  + p)  = 0 (53) 


In  \ppendix  VII  we  show  that  the  most  general  solution  to  153)  is  given 
by 


' V l 

fly)  = l fn“  - cos  ";)/)  2 - sin^j  Jl:^ 


s i n 


2lf> 


Ott 


cos 


-V  X 


(54) 


when 


re  1 is  an  a rb it rarv  even  periodic  function  of  vb  and  F is 
c t Ott 

an  arbitrary  odd  periodic  function  of  i p,  both  of  period  tt.  Thus  any 

solution  for  t lie  diffracted  field  f(tj0  must  satisfy  154). 


5.  4.2.4  Solution  using  asymmetric  contour  y . 

In  the  previous  section  we  raised  the  possibility  of  starting  with 
a contour  other  than  y and  obtaining  a solution.  Such  contours  do  exist 


and  the  asymmetric  contour 

> = y -y  , as  shown  in  Fig. 3.9, 

is  just  one 

of  these.  If  we  write 

E(p,<|>)  = 

Tvr  J 

[ iicOc'-ipcos(a-^da 

(55a ) 

and 

Y+4> 

I j (P.VM 


r 


i 


- i oncosis -ij>)  . 

' U 


( S5b) 


i 
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with  r being  the  mapping  of  7 under  the  transform.!  son  (ou.b  , we 
obtain  the  following  set  of  functional  equations  for  s and  ; 

3 IT  - Z-rr 

*sU’t+  — ) = - s(-ot+  ) 


s,  l.i(ol-  j|  = - sL  l - - (.Ot)  -n/J  J 

* *1-*)  = na)  ^[ctoO]  + s;  |- 

" st-a)  = st  |r.(oO  ] - Sj  [-z  (a)  | 


.(.aoe  1 


1 he  m°st  Ren era l solution  to  (5oa-5od)  is  giver 


~*w  = + - f)  ♦ '"offCosi—  - *, 


where  d,..le)  is  any  function  satisfying  (S7bl  with  a period  Jit  in 

'e:  t'‘1,  l07t*-a1  are  anv  evon  -l,ui  odd  functions  of  a with  period 
••  However , such  a solution,  apart  from  having  exploding  branch  cut 
integrals,  gives  rise  to  incoming  diffracted  waves.  To  see  this  we 
!l',‘  1'Kt‘  "*  ‘ lts  equivalent  sum  of  the  contours  p^,  p ( p g . ,inj  _|> 

• shown  in  l igo.IP.  In  doing  so  we  pick  up  the  residues  of  poles  situated 
between  the  strip  .n^<Rea<ir^,  and  we  may  write  the  integral  over  7 as 
the  following  sum. 

_ I'  -UHOsia-'t‘,'da  = :,tj  £ Rosilluos  + | ^aU,--'oeosu-sM  ^ 

VB„ 


* I ♦ : j 5<c,)c'jpcosl01''<'1 

n *n.  1 


id  (581 


lunation  i.-')  contains  branch  nil  interrals  i'vit  K,  and  B . 1'he 

i ii  t im;  ra  I s ovor  the  upper  half  of  B and  t ho  lower  half  of  h diverge 

as  p becomes  laryto.  Pue  to  the  different  forms  of  s(o)  and  s(x) 

it  i s not  possible  to  mutual  1'  eaneel  the  divergent  hraneh  out  integrals. 

further  the  solution  from  si  O contains  an  integral  over  the  path  1* 

1 0 

which,  when  evaluated  usinc.  the  steepest  descent  method,  can  he  shown 
to  have  the  form 


. - 10,  — 
( ■ e v i 


which  represents  an  incomini;  diffracted  wave,  and  can  form  no  part  of 
a plnsically  meaningful  solution.  Thus  the  use  of  the  asymmetric  contour 
) introduces  additional  d i f t'icul  t i es . In  the  next  two  sections  we 
examine  alternative  contours  of  two  other  types  hut  discover  that  am 


solution  obtainable  hv  usim;  these  modified  contours  has  identical  form 
to  either  of  the  solutions  we  discussed  before. 


s . I . . S Shifted  symmetric  contours. 


Y h/-nut)  * i>  dnnt) 


IVe  note  that  any  solut  ion  obtained  by  integral  i m; , s(p)  of  the  general 
term  i ISa ) and  the  co  r respoml  i m;  ('unction  s ( (a)),  over  any  arbitrate 


contour  x where  \ is  svmmetric  about  the  origin,  will  xatisfx 

a ’a 

Maxwell's  eipiatious  and  the  houtnlarv  conditions.  We  could  choose  y to 

be  ■y  where  \ is  shown  in  liro.ll  for  m 1.  let  us  examine  the 
m m 


so  1 ut i on  I 


lp..{>)  over  such  a contour. 


= slm\a)e'jpcosta-^ 

. 


da 


Y.' 


= { S(m\a+:nm)e--ipC0s^-^da  + 


slu)  (u-2miT)e'jpC0S 1 


Y ' 


C«('l 


where 


*‘"’w  • * -r-  (tt  - t)  * if  - Y)  U. 

and  primes  on  Ym.  Y+  and  y_  indicate  a shi ft  by  j . 


I ) 


By  noting  the  periodicity  of  the  functions  ) i'  (ni)  p (m1 

-'it  ’ eir  * o it 

t and  making  use  ot  the  decomposition  relation. 


fl  + I f’  - > 


(fi+f2)  + T I 


) (('-’1 


Vrj 


V^J  are  arbitrary  functions  integrated  over  arbitrary  contours 


-Trjfl'"\p,4>)  = [ l d (Cl) 


Y+vH 


1 1 ’ I j we  can  write  1:tni\p.4>)  in  the  following  form 

♦ cos Jm t [peJm>cos  (■£  - Ar) 

)]J  X e-JPCo»(a.*lJa 

*J  •">*'  - f)  j'-'V 


5.  Til 

JJj 

X 0-jP^s(a-4>)da 


where  y and  > 


are  symmetric  and  asymmetric 


is  of  the* 


contours  that  wo  discussed  earlier.  We  note  that  t:  111  (p,4>) 
f o rm 


. _ f -.ipcos (ot-ip)  - -jpcos(a-ii))  , 

'.'PI  stale  da  *■  stale  v v'da  (ivl) 


wh i 1 c*  s and  s are  of  the  general  form  (48a)  and  (5~a)  respectively, 
llms  an\  solution  over  may  be  represented  as  a combination  of  solu- 

tions i>ver  i ’ and  ) ’ . But  we  have  already  noted  that  any  combination 
ot  solutions  over  )'  and  ) 1 will  have  diverging  branch  cut  integrals 
and  or  incoming  dit  traded  waves.  Hence  any  combination  of  solutions  over 
yjn<  Y ' • ) 1 will  have  the  same  difficulties  and  fail  to  give  a phvsicallv 
meaningful  solut ion. 


l.-'.u  Shifted  asymmetric  contours, 

) m = (Y  + + dm  it  1 - tY_  - -’hit) 

In  a manner  similar  to  the  previous  section,  anv  solution,  over 


m ’m  T ’ ■ 

ot  solutions  over  the  contours  y’  and  y'.  Hence  a solutii 
a contour  does  not  possess  any  additional  advantage. 


ovo: 

v an 

! as  ; 

i sum 

over 

such 

>.  I.J.7  Contours  with  the  end  points  separated  by  Jiii't  with  m>  1 

All  such  contours  may  be  decomposed  as  a sum  of  the  form  y + y 

m 

°r  Y + \n  depending  upon  their  symmetric  or  asymmetric  nature  as  shown 
in  I ig.-xld  for  the  case  of  a symmetric  contour  with  m 2 . 
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Thus  we  have  exhausted  all  possible  symmetric  atul  asymmetric  contours 
which  do  not  cross  t he  real  axis  and  whose  end  points  at  infinity  start 
in  one  shaded  region  and  end  in  another  shaded  region.  The  end  points 
of  any  moving  contour  must  lie  in  a shaded  region  at  infinity  so  that  the 
Sommer fold  integral  does  not  alter  as  the  contour  is  shifted  with 
Further  any  contour  whose  end  points  at  infinity  are  in  the  same  shaded 
region  are  trivial  since  any  integral  over  them  would  identically  vanish 
unless  they  enclosed  singularities  by  crossing  the  real  axis.  \ny 
contour  that  crosses  the  real  axis  is  not  useful  in  our  method  since  such 
a contour  would  cut  a branch  cut  as  it  is  shifted  with  if- 

We  consider  two  other  possibilities.  One  is  to  extend  the  branch 
cuts  up  and  down  wards  to  infinity  rather  than  joining  them  pairwise. 

This  is  definitely  not  workable  for  our  present  geometry  since  such  a 
definition  of  branch  cuts  would  restrict  the  movement  of  the  contour  to 
less  than  tt  whereas  we  require  it  to  be  moveable  by  .er/d  without 
cutting  branch  cuts.  However  we  explore  this  possibility  in  the  context 
of  a different  problem  where  a movement  of  n is  sufficient.  The  second 
possibility  is  using  finite  and  fixed  contours.  We  discuss  both  these 
approaches  in  the  next  section,  and  conclude  that  both  of  the  approaches 
fail  to  give  the  desired  solution. 


.>-46 


ECP.W  • X / «<♦<? 

“ J y+$ 


(foSa) 


'!»•«  HI  J.  V >• 


ncos  i .'  -4>) 


— — < 4>  < 0 (i>5b) 


whore  the  variables  \,^  and  the  contours  y and  T have  the  same  meaning 
as  before. 

The  boundary  conditions  are  given  by 


H.P,7t/2)  = 0 


E(p,0)  = I^lp.O) 


3K(p,4Q  _ 3hltP,<{>) 


(6t>c ) 


and  (boa) - (Obd)  lead  to  the  following  set  of  functional  equations  for 


s(a)  and  Sj  (.£(a)) : 


s(a)  = sl-a+Ti) 


SjtCla))  = -Sj  (w  (-oi-iT )) 


(b7a) 


lb7b) 


s (a)  - s (-a)  = T(ot)  Is^cia))  - s^^i-ot))] 


s(a)  + s(-a)  = Sj  (r, (a) ) + Sj  (C(-a)) 


Proceeding  in  a manner  similar  to  that  shown  in  \ppendix  V we  obtain 


the  following  general  solution  to  slot). 


at , t , 


at  , t , 


s(a)  = 1’  cos l — - - — f-)  + P sin( — — - 

OTT  IT  riiT  v Tr 


where* 


cost  , = -j ; 0 < t , < IT 

- l + x(a)  - - - 


and  P and  P have  the  same  meaning  as  before. 
eTT  O IT 

If  we  proceed  by  using  the  t-transform,  as  shown  in  Appendix  VIII, 
we  obtain  the  following  particular  solutions  to  s(a)  and  s^(<;(a)). 


, ...  V„  - a {(i  + a - ir 

2 sin{ (it  - ti)}  sin{— (tt  - t •>) } 

s(a)  - -4 * 

smt,  smOJ>  -a)  sinOpg  + a) 


*0  + a 4’0  - a - tt 

sin{ — - — (tt  - t , ) } sin{ (ir-ti)) 


sin  (t(i  0 + a) 


sin(V0  - a) 


(C(a)}  = 


+ a - a - * 

, , _ . T s till ; — (tt  - t->)}  s in{ (it  - t ■>! } 

1 + T (q)  TT  IT 

2t(a)  sin(i(^  + a)  sin (tf1^  - ot)  1,111  . 


s(a)  may  be  put  into  t lie  following  form  so  as  to  conform  with  (bSa) 


s (Cl)  = 


SitT-'lTpCOS  ( - 


t 1 at , t , 

-f-)C0S(-y-  - y-) 


<|T0t . t,  at , t,  J 

• 2sin2asin( — - - — - ir)  — ; s— 

TT  - TT  - COS-.  , - COS 2a 

0 

t|»0t  J t,  at,  t , 

- ’si'H — — - -y-)cos(-y  . _^) 
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Similarly  s^C^fa))  may  be  transformed  into  the  following  form. 


•.«<»” -TisF 


fsin2t|)0cos  (- 


^ •>  ^ i t 


- . _ 'I'-.t,  t,  at., 

— )sin(— — + -f-)  - sindasinf— — --^ft)cos( — - + 

JZ ~ 7T  J 7T 


cos2^0  - cos2a 


^0^2  at-> 


si,H— ^ - -j)sin(-jpt  + -^iVsint,, 


(70b) 


1- roin  (69b)  and  (/Ob)  we  obtain  the  following  two  equivalent  expres- 
sions for  s^a): 


s: (a)  = 


f..  .Vd  t2,  . “t^  t'>  ~ ^0t->  t'1  at->  t-> 

siii-Vqcos(  - 2“)sin(-^—  + -*f-)  - sin2asin(—^-- -^-)cos( — - + ^) 


2 

(cos  \p  - cos“a) 


V?  t,  at.  t, 

2sin(— f)sin(—  + -f-)} 


1 + x(a) 

— sint, 

-x(a) 


(71a) 


where 


a 


-1 

cos 


(ncosa) 


(71b) 


t,  = t , (a) 


(71c) 


Sj(a)  = sint. 


sin{ 


Va 


-("  - 1->)} 


sin{- 


K " 


-(it  - tj) } 


sin  (t^0  + a) 


s in  (i(,0  - a) 


1 + t (a) 
2t  (a) 


(72) 
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It  should  be  noted  that  a = cos  * i is  not  a pole  of  s ^ (a)  since 
sint,  also  vanishes  at  this  point.  We  write  H and  as  sums  of  three 

terms  corresponding  to  pole  contribution,  diffraction  and  branch  cut 
contribution  as  below. 


+ Ejip.^)  + Eb(p,4i)  (,73a! 

lijte.P)  = F.  ! (p,4>)  + EdlCP,<fO  a fbl(.p,$)  C3b) 


3.  5.1.1  Pole  contribution 

Io  find  h (p,4>)  we  consider  the  following  two  cases: 


i)  ° < Ss  < ipQ  < n/2 


In  this  case  the  poles  at  -n  + i/^,  -(ft  and  n - i/ . contribute  to 

the  solution  with  the  residues  cost  ,0, -cost ^ 1.  and  -l  respectively. 


when 


"20  = VV 


(74) 


ii)  0 < iji  < 4>  < n/2 


In  this  case  the  poles  at  tt -iji  , 4'  , and  rr-41 , cont  ribute  to  the 

u 0 o 0 

solution  with  the  residues  cost  -cost 1 , and  -1  respectively. 


"20 

Thus  in  both  cases  1!  is  given  by 

P 


E (p,$)  = e-jpcosW-ip0)  _ e+jpcos(W»0) 


1 ~ 1 f -jpcos(4>+4)0i 
l + tU'0) 


- e+jpc°s^-^}.  (i  < 4>  < jt/2 


(751 


3-50 


Similarly  we  consider  the  pole  contribution  from  The  only 

poles  of  s^(oi)  that  contribute  to  the  solution  are  at  a = and  tt  - 

both  of  which  have  a residue  of  2sin.  /[/  n“  - eos'ij),,  + sini!>,]  and  E , 

o U u pi 

is  given  by 


Epl(p’*) 


[e_iPncos  1 41  - <*•'  o ) + e+ j pncos  -J  ^6a-) 


n"  - cos"p0  + sinp() 


V„  = cos  

0 n 


o < vQ  < it/: 


3.  5.  1.2  Diffracted  field 

Excluding  the  shadow  boundary  regions,  where  ij  - ^ or  t ' tyg,  we 
obtain  the  following  asymptotic  expressions  for  the  diffracted  fields 
E , and  E . , : 

U ul 


Ej(p,i,M  = -{s(ir+4>)  - se-TT+p)}  e 


j (p  + T1 


/vdtrp 


H'o  - 4»  *o  - * - * 

r o _ . . r u 


- j (P  + j)  (cos{ — - — (tt  - t )}cost  N + cos{- 


(^-t  )} 


sin  W'q  - ^ ) 


m a * - tt 


<t'0  * * 


cost (tt  - t. . 1 }cost  + cost — - — (tt  -tit 

7T  4*  P TT  (p 


sinfv'o  + ^ 


(T-al 
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■ i tT. 

J (np+T)  

F-j  (P.'t)  - ~{s^l<j)+TT)  - s j (cJ>-tt)  }e  //Jimp 


dsin“t 


1-cost 

4> 


>s{ 


ip0-<p-n 


%+<P 

^’V}  cos{—  Or-t  )} 


sin(V^ 


sin  (<|>0+<{>) 


;*?"V  ±tt-^ 


i (np  + 

4 , r; 

e //2TTnp 


(77b) 


where 


U = t2(iP) 


t = t,(<|>)  = t,(4>) 
9 


cj>  = cos  *(ncos0) 


(77c) 


(77d) 


(77e) 


If  we  let  the  refractive  index  n to  approach  unity (7'a)477b) take  the 
identical  form  given  by 


F .(p.4>)  I = E,.  CP , 4>) 


1 


1 


.1  (p  + J 


n=  1 


dl 


n=l 


’Pq-'P 

cos( — — ) s in  ( — xr~). 


o < 4)0  < <j>  / ± \i>0 


//:np 


(78) 


5. 5. 1.3  Branch  cut  Integrals 


The  branch  cut  contribution  to  the  fields  in  the  region 

0 < < rr/d  is  given  by  (Figs  .3.14a,  14b) , 


gure  1 I a . closing  of  t ho  contour  y by  means  of  the  contours  P and  P 


I m a 


>;ii  «'o  ,>  .lib.  Branch  cut  contours  H , H , It  , .itul  U in  the  plane  of 

0 < 0 :i  i ll  - 1 

the  complex  variable  o. 

Ihe  integrals  around  B and  It  are  well  behaved,  but  the  integral 

around  H.  and  It  are  both  divergent  at  infinity  as  discussed  before. 

Ot  a 

In  sec  t ion.t.1 . we  considered  all  possible  solutions  over  infinite 
contours  and  tailed  tit  obtain  a solution  free  of  diverging  waves  for  i-v 
problem  of  rectangular  wedge  on  a semi  infinite  plate.  In  the  next 
section  we  explore  two  other  methods  of  modi  f’v  i nr  .avadskii’s  solution 
1 1>  the  present  problem. 

5...’  Attempts  to  correct  i'avadsk  i i 1 s solution 

In  the  first  part  of  this  section  we  aim  to  find  a solution  (l.ljl 
obtained  by  integrating  over  fixed  finite  contours,  which  when  added  to 
.avadskii’s  dilution  would  exactly  cancel  the  first  two  branch  cut 
integrals  in  l''M.  In  the  second  part  of  this  section  we  explore  the 
possibility  of  obtaining  a solution  by  re  defining  the  branch  cuts 
vertically  to  infinite. 


• 1 Seeomla  ry  solut  ion  usim;  fixed  finite  contours  around  the 


branch  cut s 

Ke  define  the  fields  f a ml  !•'  as  follows. 
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where  t lie  contour  - of  integration  are  shown  in  lira  11c. 


iSi'al 


icucos  ( 


- i pnoos  ( 


(SOM 


I'i  gureS.l  lc  . Branch  cut  contours  B,  , B,  , B , ami  B in  the  plane 

il+  i)-  -a  + -a-  1 

of  the  complex  variable 

It  should  be  noted  that  all  the  intervals  except  those  over  R,  and 

1 0+ 

B are  convergent.  The  intervals  over  B.  and  R appear  in  l:.(p,A)  in 
such  a manner  that  the  sum  T + li  will  be  free  of  any  divergent  waves. 
Since  all  the  branch  cuts  in  '.-plane  are  along  the  real  axis  1 ^ and  1:. 
will  be  free  of  any  divergent  waves.  Now  we  examine  the  possibility  of 
finding  the  functions  l:,  !•',  t » , G,  h and  h such  that  1:  and  lij  satisfy  the 
following  boundary  conditions. 


h(p,n/2)  = 0 (Sla) 

f:(p,0)  = H. (p,0)  (Sib) 

1!’  tp,0)  = T.{  CP,0)  (Sic) 

t:.'  (P,-Tl/-)  = (Sid) 


By  a suitable  change  of  variable  we  can  write  li,  1' ' . lij,  and  I j as 


f o 1 1 ow  s . 
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-TTjn(p,<t>) 


= - J'slak-J,Y0Sl‘-;’\l.u  Js(-a+7r)ejpc0sCa+<j)5da 


0+ 


0+ 


f [h  C-ot)  + lu-'-'OJe'^00^1  ‘*',Mda  + J [hta+tr)  - h (a+7r)  Jej  l'x‘  '^da 


0+ 


0 + 


(82a) 


(-—  )H’Cp,4>)  = - vf  sin  (a-if>)s  (a)c  J- L°-s  <-lX  ■ \ia  + J*  sin  (a+<j>)  s (,-a+u  )<?-'  ‘ " °s  U v ’da 


0+ 


’0+ 


♦ Is, 


sin (a+ifO  Lh(-ix)  +h(-a)] 


e"  j Pcos  (a+lt>)  ja  . f sin(a.u  [h(a+iT)  - h(a+1T)lejpCOSl0l'^da 


0+ 


(82b) 


0+ 


0+ 


0+ 


(82c) 


w 
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nsin (;-4>)  I l'(0  + 1 


t;)lo-jp„c°s(C-«d. 


- J nsin^+4))  - Ft-^-TT)  ]cJpncos  *•  + '1d: 


nsin(.£+i{>)  [cl(-rj  + G(.-?)  ] 


-jpncos 


0+ 


J'  nsin  [G(c-ir)  - G(;-n)  ]o*’ tMKOs 


0+ 


(.8:di 


27Tj  I- 1 (p, 0)  = J (f  Cot)  + f (a))e  J^c0-S'xja_  J*  [ f (-a-n)  - f (,-a-n)  ]e-’ °C0's'\k 


0+ 


X [g C-ot)  + g(-a)]e  jtK0-saJa+  J*  [gia-ff)  - g(a-7r)  ]e-'p''0s‘\lot 

B 


0+ 


0+ 


(.88 1 


where 


f(a)  = I:(Ox(a) 

(84a) 

7(a)  = F(Ox(a) 

(84b) 

g(a)  = G(r,)t(a) 

(S4c ) 

g (a)  = (l(r;)T(a) 

(84d) 
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Similarly  we  may  write  l-:'  (p,0)  as 


0+ 


- 1 


si  not.  . . jp  cosot , 

— -ylt(-ot-iT)  - t(-a-iT)  )e-  da 


0+ 


+ I TTTT[«(-a)  +S(-a)]e'jpC0Sada 

V 


I 


sina,  , v — , , . jp  cosa , 

— ylyCa-iT)  - g(a-iT)  )eJK  da 


0+ 


1 


(85) 


The  boundary  conditions  (,31a)  and  Sid)  lead  to  the  following 
functional  equations. 


h(-a)  = h(a+iT) 

(Son) 

h(-a)  = -h(a+it) 

(Sob) 

f(a)  = -f(-a-ir) 

(86c) 

f(a)  = +f(-a-Tr) 

(Sod) 

g Cot)  = -g(-a-u) 

(Soe) 

g(a)  = g(-a-TT) 

(Sof) 
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The  boundary  conditions  (Sib)  and  Sic)  at  0=0  lead  to  the  follow- 
ing functional  equations. 


-s(a)  - h(-a)  - h(-a)  = f (a)  + f(a)  - g(-a)  - g(-a)  (87a) 

s (a)  + h(a+n)  - h(a+rr)  = -f(-a-n)  + f(-a-tr)  + g(a-tT)  - g (a-rr)  (87b) 

-s(a)  + h(-a)  + h(-a)  = — p-r-[f(a)  + f (a)  + g(-a)  + g(-a)  ] (87cl 

T fot  9 

s(a)  - h(a+rr)  + h(a+u)  = -----j-[-f  (-a-it)  + F(-a-iT)  - g(a-ir)  + g(a-iT)]  (87d) 


Combining  (87a)-(87d)  with  (86a)-(86f)  we  obtain 


-s(a)  - h(ot+ir)  + h(a+n)  = f(a)  + f(a)  + g(a-ir)  - g(a-tv)  (88a) 

s(a)  + h (ot+  n ) - h (a+ir)  = f(a)  + f (a)  + g(a-ir)  - g(ot-ir)  (88b) 

-s  (a)  + h(a+iT)  - TT(ot+'ii)  = [f  (a)  + f (a)  - g (a-n)  + g (a- it)  ] /t  (a)  (88c) 

s(a)  - h(a+fr)  + h (ot+rr)  = [f(a)  + f(a)  - g(a-iT)  + g(a-ir)]/T(a)  (SSd) 


The  structure  of  equations  (88a)- (88d)  is  such  that  they  represent 
four  equations  in  three  unknowns  and  they  are  incompatible.  Thus  we 
cannot  find  f,  f,  g,  g,  h,  and  h such  that  E,  E would  meet  the  boundary 
conditions  (81a)- (SKI).  However  the  boundary  conditions  (Sib)  and  (Sic) 
may  also  be  satisfied  by  requiring  that 
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j f siiials(a)  - h(a+ir)  + lUat+tT)  ] sin  (pcosa)da 


J [t'(a)  <-  f(ot)  - g(a-ir)  + g (ot-n ) | -^ycos(pcos  d (91b) 

Bo+ 

Because  of  the  nature  of  the  path  of  Integration  it  is  difficult  even  to 
answer  the  question  of  existence  of  a solution  to  these  equations.  thus 
we  seem  to  have  come  to  a dead  end.  In  the  next  section  we  investigate 
the  possibility  of  obtaining  a solution  by  extending  the  branch  cuts 
vertically  to  infinity. 

3.  5.2.2  Formulation  with  infinite  branch  cuts 

In  this  section  we  formulate  the  problem  by  extending  the  branch  cuts 
vertically  to  infinity  as  shown  in  Fig. 3. 15.  We  define  F.  and  F,  as  tallows 

2irjF:(p,4>)  = /s(a)e-jpC°^a^Mda  0 <■  i}>  < tt/2  (.‘d-'a) 

Y+4> 

2 n i F. . (p,<j>)  = J s U)e  JpnC°sU'lMdr.  0 > f > -ir/2  (92b) 

l'+4> 


who  re 


-1 

cos 


COSO 


) 


I'd  2 cl 


and  1’  is  the  mapping  of  y into  2,-plane.  With  such  a definition  ot  branch 
cuts  ;,(oi)  and  Hot)  satisfy  the  following  relations 
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r 


5(a)  = a-oO 

(93a) 

5(a+tr)  = -5(a)  + tt 

(93b) 

i(a)  = -t(-a) 

(93c) 

T(a+TT)  = -i(a) 

(93d) 

As  shown  in  Appendix  IX  the  boundary  conditions  of  the  problem  lead 
to  the  following  functional  equations  in  s and  s^ 


s (-a)  = s (a+u) 

(94a) 

SjC-5)  = -s1(-5(a+ir)) 

(94b) 

t(a)  [s(a)  * s(-a)]  = s1(5)-s1(-5) 

(94c) 

s(a)  + s(-a)  = sx(5)  + s^-c) 

(94d) 

If  we  try  to  solve  these  equations,  as  we  have  done  in  Appendix  X, 
we  obtain  expressions  for  s^(5)  and  s^(-C)  which  are  not  compatible  with 
each  other. 


r 
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5.0  CONCLUSIONS 

In  t n report  we  have  made  a thorough  investigation  of  Zavadskii's 
method  to  tr>  to  obtain  a solution  to  electromagnetic  diffraction  problems 
involving  a rectangular  dielectric  wedge  (0  ^ <j)  -tt/Z)  and 
i)  infinite  metal  plate  along  4>  = ± tt/2 
i i 1 Semi-infinite  metal  plate  along  = - r/ 2 
i i i 1 Perfect  magnetic  conductor  along  if>  = - n/Z  and  a semi- infinite 
metal  plate  along  4*  = ^/-. 

In  all  the  cases  that  we  considered  we  found  that  Zavadskii's  method, 
as  it  is,  gives  a solution  involving  branch  cut  integrals  that  grow 
exponentially  in  the  far  field  thus  violating  the  radiation  condition, 
for  the  trivial  case  involving  an  infinite  metal  plate  we  have  found  a 
simple  way  of  modifying  the  solution  so  as  to  conform  with  the  known 
exact  solution.  We  have  made  several  attempts  to  modify  Zavadskii's 
method  to  obtain  a solution  satisfying  the  radiation  condition  for  the 
cases  i i 1 and  i i i 1 above  but  none  of  our  attempts  proved  to  be  successful. 
However  one  of  the  methods  involving  a secondary  solution  with  branch 
cut  integrals  alone  lead  to  two  integral  equations  whose  solution  seems 
to  be  either  very  difficult  or  impossible.  In  the  process  of  these 
attempts  we  discovered  the  solution  to  the  quasi-trivial  problem  of 
illuminating  a rectangular  dielectric  wedge  resting  on  a semi- infinite 
metal  plate  with  plane  waves  such  that  there  is  no  net  diffracted  wave 
from  the  edge. 
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over 
in  teg 


where 


In  this  appendix  we  obtain  an  expression  1^  for  the  integrals 
the  branch  cuts  + and  ' =Bn  + and  s*'ow  that  these 

rais  grow  exponentially  as  p-*»t 


lb  = f s(a)c-jPcos^-^  da 


s(«)c--ipcos(a-‘1’1  dc 


...  -jpcosiu-<i>)  , t 
s(a)e  v * da  + 

B 


s (a)  e”'ipcos  (a'^^  da 


TT- 


IA1 . 1 ) 


'\)+  "()-•  btt  + ’ ;l,ul  B are  shown  in  Fig.  Al.l. 
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By  a suitable  change  of  variable  we  may  write  1 ^ as 


lb=  j s(a)e 

0+ 


■jpcos  («-<(>) 


da  - s(-a)e-jpcos(a+^ 


jpcos(a-ct>)  . f , . jpcos(a+<i>) 

sla+TT)eJ  v T da  - s(-a+ir)e  v ^ ’ 


da  (A1.2) 


Noting  that  s(-a)  = s (a+rr)  we  may  write  1^  as 


Ib  = | s (a)  [e-jP«s(a-t)  _ ojpcos(a+#)J  Ja 


j"e-jpcos(a+<j))  _ e_ipcos  (a-4>)J 


I , . -jhsina 
-j  s(a)e  sin(acos 


(,acosa)da  + 2j  s(-a)e-'  1s*nct  sin(.acosa)da 


where 


(A 1.3a) 


a = pcos4> 


b = psin4> 


s (a)  = 


,.  . cosip., 

M-sina  0 


sina-  sintg()  M+sina  sina+sin^ 


(A1 ,3b) 
(A  1 . 3c ) 
(A1.3d) 


2 , H 

M = (n“-cos“a) 


(A1.3e) 


and  the  sign  of  M must  be  chosen  as  shown  in  Fig.  A 1 . 1 . Noting  that 
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only  the  second  term  in  (A1.3d)  will  contribute  to  the  branch  cut 
integral  and  choosing  the  proper  sign  for  M we  may  write  1^  as 
the  following  sum  of  two  line  integrals. 


J3, 


lb  =-8jcosil>0 


| M | s i not 1 

n-_!  sina+  sin^0 


- ibsina  . . . 

e • sin(acosa)da 


j Pi- 


Sjcos(|/ 


0 


s i na  1 


2 . sina-sinty- 
n - 1 0 


ibsina  . .. 

eJ  sm(acosa)da 


*0 

= Sjcos^ 

0 


1 Nl|  sh3  1 
n“-l  isW+sin% 


ebsh^sin(ach3)(i3 


- 8jcos4>0 


0 


0 


| M|  shll 1 

n2-l  lshe-sin1'o 


- ItshS 


sin  (.achpldB 


(A1.4a) 


where  ch3„  = n (A1.4b) 

Since  a and  b are  positive  when  0 < if)  < tt/2,  it  is  immediately 
clear  that  the  first  integral  ( A 1 .4a)  grows  exponentially  as  p-*00  while 
the  second  integral  converges  for  all  p.  lie  note  that  the  diverging  integral 
corresponds  to  the  unshaded  branch  cuts  B(1  + and  . 


I 


a 
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Al’lT.Nl'  1 \ 11 


In  this  appendix  wo  obtain  s(a)  as  given  by  (33a)  by  inverting 


t(.x,31  which  is  given  by  (32) 


. , , TT  ,, 

eh(^i'l-r  (a)  snivel 

tut  ,31  = 2sh  (,  1 ,-3n/2)3  • -p; - — n — ,-r~ — rr “TT ; — r (.\2.11 

v 0 | 1 t-v  (ct)  ] shditis*  [ 1 -T(ot)  1 simp 


s (n ) - t (oi , 3 ) e d3 


(A2.2) 


s(o)  may  also  be  written  as  the  sum  of  the  following  integrals. 


(o)  1 1 (a.31  +t (a, -31  lclux3J3  - 1 1 (oi,3)-t  (oi,-3)  |shot3d3  (A2.31 


After  substituting  for  t(a,3i  and  t(Gt,-0)  we  obtain 


s!i|l^0-3n/:ib|  |ch(^)ehoc3+T(a)sh~3sha3] 

" “ [1+T(al]  [ehno+A) shrr3  1,0 


(A2 . 4a 1 


a - 1 l~U«)  _ 

2 l*T(oi) 


:ost ; 0 < t < it 


(A2.4b) 


Noting  that 


|eh  '.v'Olvv*  i (oil  sh^bshob  | 1 1 ♦ r(a')  Jchls  ♦ixl3*  1 1 - 1 (nl  |ch(  ! -o03 


(A2.S1 
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WO 


write  s(oi)  as 


S(o) 


shlu,l)--'n  J1  (5 1 (eh t ^ nxH'  + J.VIu1’  -ah' I 


(clrv*  \)  slrv 


<33 


I | sh((ip0-tr+a)8]-*sh[  (ij)Q-2n  *J  6]  ♦ 2A{sh  f C^0“W-a)3]  *sh  [ (^0-2tr*ot)<a J > 

J j (chno‘  Alsh’tH 


. do 


(A 2.M 


winch  can  he  immediately  put  into  the  form  (33a). 


\rn\ni\  m 


In  this  appendix  we  transform  s(o)  as  given  by  (33a)  into  the 


f o rtn  ( 3o  i . 


s(a)  - 1 1 1*  +a-n)  ♦ ,Uv  ,-ci-Jn)  ♦ \I  U1  .-o-nl  + VI  (V.,+o-.2ii) 


l A3. la) 


HtM  • - - Jjsi.r  ^rll  (l-A)sin*  ? 

(1-A  )sin0  1 " ' 


( 1 - A“ ) s i nO 


jl  • \ 1-  \ 

! ■ ' 


cos  (0-6t  it) 


(A3. lb) 


A ■ -COS  t 


(A3. 1c) 


Ho  note  that  1(0)  has  a periodic  and  a nonperiodic  component . IVnoting 
s,„  as  the  periodic  part  of  s(o)  we  obtain 


i/:  -i/:  1/: 


Jn  (1-A)  Isiiu.  siny.  sin'*,.  sinv. 


. COSv'  , cosy , cosy  , COSV  . 

^ f 0 ^ * U^y  S 

(1*A)  \5  s iiu  2 S ills  . s ir.v  S1IU  . 

1 S U M > 


(A.3.:al 


where 


7S 


3-  7b 


where  we  have  made  use  of  (\5.1c)  and  decomposed  the  cosine  products. 
Upon  further  decomposing  (Ad . 1 ) we  obtain 


s(a)-s,  (it)  =- — — ^'situdv  ,-Vst/T+5t)  - sin(2a-ij' .t/u+5t) 

( 1 - A) 


s i n C - v - v jt/Tr)  + sin(2a+v^t/Tr)S  (cosda-cosde^) 


<sin(2^0-i/»  t/ir+3t/2)  eos(at/ir-3t/2) 

f 1 - A 1 


cos  (.Ja-cit/Ti+dt/dl 


siiuv  t/tT-5t/2)j>  (cosdot-c 


os. 


~ |[sin:^0cosiv0t/TT-5t/:i  -cos2t))0sin(ii'0t/-ir-dt/2)  ]cos(at/n-3t/2) 


U-A") 


[cos2acos  (at/u-3t/2)  + sin2asin(at/ir-3t  2)  ] s in ('f^t/ir-dt/d  "^/(cos-la- 


-1 


^ < s in2\|)_cos  (’Kt/ir -3t/2) cos  (at/ it- 3t / 2)  + 

U-A“Hc'os2a-cos2v0)  L 


s in2as  in  (.v  t/tr-3t/2')  s ln(ott/ir-5t/2) 


w\ 


— T-  s i n (.v^t/ii-  At  / 2)  cos  (at/tt- 3t/ 2 ) 


(A3. 51 


1-A 


1 


cos2i|<0) 


By  adding  (A3. 3)  and  (A3. 5)  we  immediately  obtain  (361. 
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Al’tTNUIX  IV 

In  this  appendix  wo  show  that  any  branch  cut  free  solution  s^al, 
to  the  functional  equations  (.lla-lld)  with  = 3^/d  and  0 = it /J, 

must  he  periodic  in  ct  with  period  Jit. 

The  functional  equations  are  as  follows: 

s(a)  = s(.-ci+3tt) 

Sj  t-M-H  1 = — -l- — l s (.at)  -s(.-cn) ) + \ [s(oO+s(-oO  ] 

-t  l ci) 

h")  = 

where  tAl.lal  and  (..VI.  Ic)  are  modified  forms  of  (llal  and  (llhl  and 


l.Vt.lb)  is  obtained  by  combining  (llcl  and  (.lid). 

let  s^(oi)  be  any  solution  of  (..\4.1a)  with  no  branch  points.  Then 

s0l-a)  = s0(oi+3n)  (A4 . Ja ) 

= s pWiO  (A4.Jb) 

Sjl  -ht))  = 1 [s  1*1 -s0px+3Tn  ] + i Is  (a)+s  Irt+Sio]  (Al.Je) 


and 

Sjl  M-cx-tr)  = y— ls0lot+4rn-sflhx+n)  | + 3I C«+4tt ) +s 0 CvX+tt')  ] l.U.Jd) 

Since  Sj  (.  ;.(a) } has  two  terms  one  involving  t (a)  and  the  other 
independent  of  t (.a)  , we  must  equate  these  terms  separately  to  the 
corresponding  terms  in  s.  (~.(-a-nl  so  that  (A4.lcl  is  satisfied.  This 
results  in  the  following  two  functional  equations  for  s(1(.ol 


t A4 . la) 
CA4 . lbl 

( A4 . 1 c 1 


which  gives 
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!'(")  - |>(a+27i)e~-*-t  = 1 

Now  wo  make  the  subst  i 1 1 1 1 ion 
P(")  = 4(ot)/(l-e'jJt) 

which  results  in 

Ml")  - i|(a+>ine'J“t  = l - 


Let 

Ml")  = 1 + k (cx ) 


(AS. 13b) 


(AS.  1 la) 


(AS.  Mb 


(AS. ISa) 


giving 

1 + k(a)  - e ■* " C 1 1 +k (h+Jtt  ) | = l-e"-'_t 


aiul  hence 

k(a)  = k(a+2ir)e"-'‘‘t 
Now  let 

k(a)  = m(o)ejat/lT 
which  gives 

ni(a)ejat/7I=  m(a*2ir)  ejat/lr 


(AS. 15b) 


(AS. 15c) 


(AS. loa) 


(AS. lob) 


whose  solution  is 

m(«)  = m(oc+2ir)  = M.  (a) 
2n 

There  t'o  re 


(AS. loc) 


s (>i)  T (a)  [MM,  (ot)  e 
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.iat/ir. 


- joit/n 


♦ ^i, (a^  1 1 *^2n  ( *)  e~Juit/n  ) jat/iT 

(l-ejJt)  1 + ":7T 


(a)  (AS. 17) 
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APPENDIX  VI 

In  this  append i x we  examine  the  branch  cut  integrals  (35b)  and  show 
that,  for  the  contribution  from  unshaded  region  to  he  zero,  the  functions 
’’err  and  p07i  "Ulst  identically  vanish,  from  (35b)  we  have 


-Xilyp.M  = j s(a)e*jpcosla“l'\ll. 


(AO. la) 


where 


I)  = 


r 15(,  f 0 < <P  < II 


8,j  + 8 , n < 1 1>  < 3 n/ . 


(All.  lb) 


and 


Bo  = 

B , + 

0 + 

V 

B = 

IT 

B v 
n + 

B 

rr- 

IK  = 
Jtt 

B2n  + 

+ Sr- 

(Ad. lc 

As  we  note  from  Pig.  Ad. I,  B(H>  are  in  unshaded  region  for  the 

range  0 < <)>  < tt . In  this  range  of  4>  we  may  write 


-"jiyp.W  = l'b  (p.ifr)  + Eb  (p,4.) 


(Ad. Ja) 


where 


"lil 


(p.«M  = | s (a,  oJ  p™*  da 


(Ad.Jbl 


B + B 
D<-  il- 


and 


w 


I'ieuro  \o . 1 Branch  cut  contours  I!  .It  , It  . It  |t  .i  mi  It 

o t o - :i  *■  " - J;i  * .in 

Hie  most  jtotUMM  1 solut  ion  to  sto)  is  j'.ivoti  by  ( ISI 


ot  .tt , ,,  , nt  ot . 

s l A 1 J.to)  * I .10)  cost  - , ' f I',,  lolsnu  >'  t 

where  J to)  is  .ut'  periodic  function.  of  pertoil  s.itisfvmr,  th 


re  lot  iv'n 


J ...  (a  * v-  > d (-a*S*  i 


l AH  . .'ll ) 


and  I'  and  P ...  are  arbitrary  even  and  odd  periodic  functions, 
Ot>  t > a * 

respectively,  of  a of  period  it.  t is  defined  by 


cos  t = 


ll-T(.al  1/ll  + ti.ci)  | 


(An.  5c) 


We  may  choose  J,  lot)  to  be  free  of  branch  points,  so  that  it  does 
not  contribute  to  H^IP,*?).  Now  let  us  examine  the  branch  cut  integrals 
in  unshaded  region  which  are  given  by  (Ao.db),  and  write 


f,.  IP.40  = 


R , + B 

0*  TT  - 


,,  01 1 
P COS  I — 
on  n 


at,  „ 

-r-)  + P.  sin 
2 Ott  v n 


at  at 


-.ipcos  la-vi) 


da 


e- .ipcos  (a-^) 


da 


1'  cos  i 


at 


r 


t1 


,+.i  pcos  (a+p) 


l An. 4) 


where  we  have  made  use  of  the  periodicity  of  P and  I'.  and  transformed 

ott  On 

the  integral  over  R to  an  integral  over  R,  . 

71  - ’ 0 + 

In  (A(i. 4)  we  do  not  see  any  possibility  of  the  mutual  cancellation 

between  the  two  integrals  because  of  the  different  nature  of  the  exponential 

terms  in  the  integrand.  The  onlv  wav  in  which  P (p,^)  could  vanish 

bl 

for  all  p and  if  is  that  the  non-exponential  parts  of  the  integrands 


be  either  identically  zero  or  be  free  of  branch  point. 


Since  t does  have 


a branch  point  wo  conclude  that  the  integrands  must  be  identically  :.ero 
which  implies  that  and  P(1J]  must  separately  vanish.  ho  arrive 

at  similar  conclusions  for  the  range  tt  < <p  s .Vi/J. 
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Al'i’i'NP  1 \ VII 


In  this  appendix  we  find  the  most  general  solution  to  the  following 


functional  equation. 


f C<j>)  - f (Tr-qj ) - RfLTT*<p)  = 0 


(A". la) 


where 


-lie  i ->  e 

R = [sinif-  (.n“-cos”ip)  ] / [s in0+ (n”-eos-i|>)  ] 


CA~ .lb) 


Substituting  for  R and  rearranging  we  obtain 


">  ■)  1* 

f(<f>)  [sin4)  + (n“-cos“4>) 2 ) 


1 1 *5 

f lsin4>  + (n“-cos“<j>)  ] 


9 *> 


[sin(ir-Hp)  + (n“-cos~(H+4>))  ] = 0 


(A-.Jl 


■>  h 

f (>?)  = g (d> ) [ in“-cos“4>)  - sin^] 


(A-. 3a) 


g(4>)  - gUM)  + g(TW)  = 0 


CAT. Sb 1 


putting  for  ^ we  obtain 


g(-4>)  - g(*+4>)  ♦ g ) = 0 


t AT* . act 


Adding  (7.3b)  and  (7.3c)  we  obtain 


5 ')  0 


SI  Cl  ♦ gt-vO  = 0 


(AT. 4 


Hence  g is  an  odd  function  ot'  o. 
Let 


gl-Jl  = J sin\i{>  fpMdX 


and  substituting  in  I'.db)  wo  got 


t\\i  | s in\vf  - s in\ ipi-ol  ♦ sin\(T+Cl]  d\  0 (A'.ol 


which  gives 


t't \1  (.  1 ‘-cos' " ) s i n\cd\  0 


llonco  cc*s\u 
and 


- V:  \ = dn  • : 


gUs1  V \ sin  diJln+1  51  + ^ B sin  d^ (,n- 1 ,'51 


which,  attor  combing  the  terms,  may  bo  rewritten  as 


gi^l  = cos  1 i'  sinde  + sin  -5?  V l>  cos dne 
•'  o " •'  5 n 


l AT . 7 ) 


l A 7 . 8 1 


(AT. 9) 


l.V.  101 


Iho  series  in  (~.10)  represent  odd  and  even  periodic  functions  of  0 
with  period  it.  So  that  g(^1  is  given  by 


w 
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APPENDIX  VIII 


I ii  this  appendix  we  obtain  a particular  solution  to  equations  (6?a)-(67d) 
using  the  t -transform. 

Defining  t(a,8),  t ^ t i,p)  to  be  the  transforms  of  s(a)  and  s ( - (a) ) 
and  continuing  the  representation  s(  <.)  = t(o,S)  beyond  the  pole  at 
ct  = ij;  , we  obtain  the  following  system  of  functional  equations  in  t 
and  t . 

^0  -Bit/'1  Rti/'i 

[ t (a  +tt/ 2 , 8 ) + e u]e  0 = It  i-a+u/2 , -3)  + e °]e  " (AS.  la) 


t (a,8)  = -t^-a.-tfje 


t(a,3)  - t (-t , -8)  =t (a) [ 1 1 (a, 3) -t  (-a, -8) J 
t (a,8)  + t(-a,-8)  = t^ci.S)  + t (-a,-3) 


(AS. lb) 
(AS.lc) 
(AS. Id) 


The  first  of  these  equations  may  be  transformed  to 


t (ci, 6)  = t(-c*+TT,-3)elf^+  de7^"  sh [3  (1t/2-4,„)  J (AS. 2a) 


Eliminating  t^  from  the  last  three  equations  we  obtain  the  following 
relationship  between  t(a,3)  and  t(-a,-8). 


t(-ot,-3)  = — — j-  t(a,8) 


(AS. 2b) 


A = cost.. 


I-TfcQ 

l+t(aj 


(AS. 2c) 


Combining  (AS. 2a)  and  (AS. 2b)  we  obtain  the  following  single  functional 
equation  for  t(ct,B)- 


-T 

t(a,B)  = ^t(a+Tr,S)  + 2e7rt3/2sh[S(^0-TT/2)]J  — 


(AS. 3) 


which  has  the  following  periodic  solution 


tfa.IS)  - sh[6(V„/2)] 


(AS . la) 


where  p («,3)  = (e_TrB-A)/ (A-e718) 
Upon  simplification  we  obtain 


(AS. 4b) 


tCa.0)  = sh[3(4,0-ff/2)]  lV^JJ?(S.^:^-.Uch(3TT/2) 


We  obtain  s(a)  by  inverting  (A8.S) 


t(a,S)e'aSdB 


(AS. 5) 


CO 

[t  (a,3)+t(a,-3)  JchagJS  - ft  (a, 3 ) -t  (a, -3 ) 1 sha$d3 

n 


(A8 .6) 


Substituting  for  t(ct,6)  we  obtain 


5 •>.. 


\m  \m\  i\ 

In  this  appendix  we  derive  the  functional  equations  (,‘d  la ) - (iMd)  . 

I he  branch  cut  structure  in  (-and  -planes  is  shown  in  1 i^.j  15.  We 
introduce  a small  amount  of  loss  factor  in  the  dielectric  constant  of  the 
medium  so  that  the  branch  points  are  slightly  displaced  from  the  real 
axis  in  the  oi-ptane.  When  the  branch  cuts  are  defined  as  in  ria.5.15, 
to)  and  Ho)  satisfy  the  follow  ini',  relations. 


• to)  t-o) 

tot'l  1 = - l (X ) Ml 

T to)  = - t t-O) 

rlo+n)  - -rto) 

We  define  the  fields  f and  I 

w rite 


l A‘>.  la) 
l Ail.  lb) 
t Ai» . 1 c ) 
t Ail  .Id) 

as  in  t'd-la)  and  ti’-'h)  , and 


-'njl.to.i,*') 


s t>x  **(> ) e 


i pcoso  . 

do 


t.V.l.J) 


for  the  field  1 , tdoM  in  the  dielectric  we  write 


-tt.'l’,  (. P . 4> ) | Sjt  a^e 


t I'llCO s 


| | s j (•'. +<f ) -s  j ttjil  |e  l,IUl  d 
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where  wo  have  made  use  of  t ho  relation  t (a)  s inot/ns  i n.. . Using  the  relations 
l'\'*  ■ - ) . (A9.3),  (A9.5),  ami  (A9.7),  and  enforcing  the  boundary  conditions 
at  vp  = O,.1  /J  we  obtain  the  desired  functional  equations.  By  equating 
the  field  l(p,4>)  to  :ero  at  if  = v/1  gives  (‘Mai,  and  bv  equating  1 (p,0) 
to  t:  (p.o)  and  I’lp.O)  to  l:|(p,l>l  we  obtain  (91cl  and  (bid) . It'  equating 

1J  IP.-tr/-)  to  zero  we  obtain 


I 
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\s(a)-s(oi+iT)  = s (.a+it) -As  (a+Jtr) 


(A  10.4a) 


where 

A = ll-Tj/d+T) 

Dividing  l. A 10. 4a)  by  /A  we  have 

/A  s(.otl  - — s(a+u)  = — s(a+Tr)  - /A  sfa+diT) 

/A  /A 

Noting  that 

A(a+a)  = 1/A 

we  have 

P(a)  -P(a+ir)  = P(.a+ir)-P(a+-tr) 

where 

P(a)  = /As(.a) 

Let 

Qlal  = P(.a) -PCa+tr) 

then 

Q Cot)  = Qla+ir)  = 0 (.a) 

where  (a)  is  an  arbitrary  function  of  a with  period  ?r. 
Substituting  (A10.76)  in  (,A10.7a)  we  have 
P(ot) -P(a+TT)  = Q (a) 

Therefore  the  most  general  solution  to  P tot)  is  given  by 

P(a)  = 1'  (a)  +a(l  (a) 
it  it 

so  that 

s (a)  = | (a ) +0(1^  (a)  ] //A 


(AID. lb) 


l A 10.5) 


(A  10 . oa) 


(A  10. bb) 


l AID. 6c) 


l. A 1 0 .7a) 


lAlO.'b) 


lAlO.Sa) 


lAlll.Sb) 


where  Fn,  (1- 


arc  arbitrary  periodic  functions  of  a.  with  period  it. 


•■>- 100 


But  s(.a)  must  also  satisfy  (Aid. la).  IV  e decompose  I and  d into 

n n 


their  even  and  odd  parts  and  write 


and 


f 1 0 I 10  * I'  lal 
rr  e o 


f>  (a)  = <;  ia)  + i;  (.o 
e o 


1 \.  id. da  ] 


(All). Ob) 


Then 


s (-a)  = S 1;  (a)  - F (.a)  - ad  la)  +,id  (a)  s A 


l MO.  10a ) 


s (a+ii)  = 


■jj^ia)  + To(al  +a[d^ia)  + d^ia)]  +rr  (d^ia)  + d^ia)]^ 


l Aid . 10b) 


Now  to  satisfy  (All),  la)  we  must  have 


2Va) 

*2a  C.e(a)+ir[t'.^(a)+G  (a)  ] = 0 

(A10. 

11) 

which  requi 

res  that 

Va) 

= 0 

(Aid. 

Ua) 

and  0 (a) 

= — F (a) 

7T  0l 

(Aid. 

12b) 

so  that 

slot)  = 

IF  (a)  + ( l-2a/ir)l:0(ci)  ]/>  \ 

(Aid. 

13) 

which  may  be  written  in  the  following  form, 


s(a)  = U’e(«)Mci-lT/’)P0la)  )//A  (.MO.  14) 

where  1'^  and  1’  are  arbitrary  even  and  odd  periodic  functions  of  a with 


period  tt  . 
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ABSTRACT 

\ detailed  study  of  Aleksandrova  and  khichnyak's  method  of  obtaining  a 
solution  to  the  problem  of  electromagnetic  plane  wave  diffraction  hv  a 
rectangular  dielectric  wedge  is  made.  The  spatial  integrals  involved  in 
deriving  the  integral  equation  by  their  method  are  systematically  carried 
out  in  full  detail  and  the  implied  convergence  criteria  are  thoroughly 
discussed.  Besides  pointing  out  some  typographical  and  expository  errors 
and  a mix  up  on  the  sign  convention  in  the  time  variation  in  their  paper, 
it  is  shown  that  Aleksandrova  and  Khichnyak  omit  some  terms  which,  if 
properly  accounted  for,  would  lead  to  an  integral  equation  that  is  not  only 
different  from  the  one  they  have  solved  but  is  also  not  amenable  to  solution 
using  presently  known  standard  techniques.  A formulation  using  a modified 
contour  of  integration  was  attempted  with  a view  to  obtaining  a singular 
integral  equation  for  the  weighting  function  that  might  be  amenable  to  solution, 
but  the  attempt  did  not  prove  to  be  successful.  It  is  concluded  that  the 
problem  of  wave  diffraction  by  a rectangular  dielectric  wedge  has  not  been 
Solved  by  the  method  under  review,  and  that  this  conclusion  must  hold  also 
for  the  later  work  on  arbitrary  wedge  angles  using  the  identical  method  - 
though  since  no  details  are  given  no  detailed  analysis  can  be  made. 
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4.1  INTRODUCTION 


l'ho  problem  of  wave  diffraction  by  a perfectly  conducting  metallic 
wedge  has  been  solved  [1,2].  The  analogous  problem  of  diffraction  by  a 
wedge  with  impedance  faces  lias  also  been  solved  [3-  10].  A generalization 
of  this  solution  to  the  case  of  dielectric  wedge,  either  free  or  resting 
on  a semi- inf inite  metal  plate,  entails  serious  mathematical  difficulties. 

Apart  from  being  a classical  boundary  value  problem  the  diffraction 
of  electromagnetic  waves  by  a dielectric  wedge  is  of  particular  interest 
in  the  theory  of  dielectric  wave  guide  matching  ] 11,1-),  radio  propagation 
over  the  earth  [13,14],  and  in  radar,  for  the  effect  of  scattering  by 
dielectric  radomos  [IS].  An  analogous  problem  to  that  of  the  dielectric 
wedge  is  encountered  in  the  field  of  acoustics  [lo]  and  in  seismo logical 
situations  involving  the  behavior  of  Rayleigh  waves  at  the  boundary  between 
the  ocean  and  the  earth  [17, IS].  It  is  not  surprising,  then,  that  consid- 
erable research  effort  has  been  directed  towards  the  problem  of  diffraction 
by  a dielectric  wedge  [19-30]. 

l'ho  problem  of  electromagnetic  diffraction  by  a dielectric  wedge,  with 
or  without  a semi- infinite  metal  plate  on  one  face,  is  a special  case  of 
scattering  by  sectoral  media.  The  specific  problem  of  the  diffraction  of  an 
E-polarized  plane  wave  by  a right  angled  dielectric  wedge,  whose  refractive 
index  is  limited  to  a certain  range  of  values,  lias  been  solved  theoretically 
by  a number  of  authors.  Radlow  [IS],  Kuo  and  I’lonus  [20]  and  Kraut  and 
Lehman  [21]  offer  solutions  to  the  problem  by  a generalization  of  the 
function  theoretic  method  of  the  Kiener-llopf  technique  [22]  from  one  to 
two  complex  variables.  However  they  do  not  simplify  the  final  results  and 


.t-;  remarked  In  Kuo  and  I'lonus  [JO]  "...the  solutions  are  too  complicated 
to  be  used  practically".  Kraut  and  Lehman  |dl]  claim  that  Radlow's  solution 
[li'l  is  incorrect.  Kuril ko  Idd.Jl]  obtains  a solution  in  the  form  of  a 
rather  complicated  system  of  Fredholm  integral  equations  which  have  to  be 
solved  by  numerical  techniques.  Fat's  [25-d~|  final  result  ends  up  as  an 
infinite  system  of  lilbert  singular  integrals  which  lie  states  are  suitable 
for  numerical  computation.  However  he  does  not  actually  obtain  any  explicit 
results  of  practical  use. 

Larj'  and  Solfre'  [I'l  have  used  an  approximate  technique  known  as  the 
Raleigh-nans-l'orn  111-11-111  approximation  |.'S|,  to  solve  the  problem  of  a 
dielectric  wedge,  whose  refractive  index  is  near  unity,  placed  on  a perfectly 
conducting  infinite  plane. 

Rawlins  formulated  the  boundary  value  problem,  of  the  diffraction 

of  an  1-  or  ll-polari '.ed  electromagnet  ic  line  source  by  an  arbitrary  angled 
dielectric  wedge,  and  obtained  a solution  in  the  form  of  a Fredholm  integral 
equation.  Using  a standard  perturbation  technique  he  obtained  a Neumann 
series  solution,  to  the  integral  equation,  which  converges  when  1 < n < >2 
where  n is  the  refractive  index  of  the  dielectric  wedge. 

\ number  of  years  ago  Zavadskii  [oO]  proposed  a method,  which,  he  claimed 
would  give  exact  analytic  solution  to  a class  of  two-dimensional  wedge 
diffraction  problems  including  the  problem  of  diffraction  by  a rectangular 
dielectric  wedge  resting  upon  a semi- infinite  perfectly  conducting  plate. 
However  the  solution  obtained  by  his  method  contains  branch  cut  integrals 
which  give  rise  to  waves  "diverging"  at  infinity  in  complete  violation  of 
the  radiation  condition.  In  an  earlier  work  |ol]  we  made  several  attempts 
to  modify  lav.ulskii's  method  so  as  to  remove  these  drawbacks  but  did  not 
meet  with  success. 
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In  section  4.2of  this  report  we  give  a brief  account  of  their  method 
and  point  out  the  discrepancies  in  their  evaluation  of  the  integr;  Is  and 
the  interpretation  of  the  residues.  In  section  4.3  we  start  with  the 
correct  expression  for  t he  incident  wave  and  systematically  reduce  the 
integrals  involved  to  the  desired  form  and  show  the  missing  terms  in  their 
integral  equation.  We  conclude  that  if  these  terms  are  properly  included 
the  resulting  equation  is  not  amenable  to  solution  using  techniques  known 
to  us.  The  finer  details  involved  in  proper  evaluation  of  the  spatial 
integrals,  which  are  not  mentioned  in  [AK-1]  are  given  in  appendices. 

In  section  4. 4we  give  a brief  account  of  our  modified  approach  to  this 
problem.  Conclusions  are  given  in  section  4.5.  In  our  analysis  we  will 
make  frequent  reference  to  their  paper  on  rectangular  dielectric  wedge 
[AK-1],  and  equation  numbers  quoted  from  it  will  be  annotated  with  the 
letters  ' AK'. 
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Maxwell’s  equations,  in  integral  form,  give  the  elect  rie  fieUt  I 
the  mar.net  ie  field  II  everywhere  \ ia  the  equat  ions 
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The  first  feature  of  the  solution  of  (d.l)  is  that  the  internal  field 
in  the  medium  (.reV)  is  directly  determined  through  the  unperturbed  field 
of  the  incident  wave.  The  scattered  field  for  (rjV)  is  determined  through 
the  known  internal  field  by  the  same  relations  (J.l). 

The  second  feature  is  that  a formal  expression  for  the  desired  field 
in  the  medium  is  the  sum  of  the  fields  of  the  unperturbed  incident  wave 
and  ot  the  waves  formed  by  the  integral  terms  [5oj.  In  evaluating  all 
integrations  this  term  should  yield  a series  of  terms  including  a wave  with 
a propagation  constant  coinciding  with  that  of  the  unperturbed  incident 
wave.  In  accordance  with  the  Oseon -Uwutd  "extinction  principle"  |o7  bp| 
this  wave  must  exactly  cancel  the  unperturbed  incident  wave.  It \ imposing 
this  condition  one  obtains  the  amplitudes  and  the  directions  of  propagation 
of  the  penetrating  plane  waves. 

T.d.d  Plane  wave  scattering  by  a rectangular  dielectric  wedge 

We  now  consider  the  specific  problem  treated  by  Aleksandrova  and 
khi.'hnyah  i.e.  the  scattering  of  a plane  electromagnetic  wave 
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by  a rectangular  dielectric  wedge  (l;ig,4.1)  with  a relative  permittivit\ 

>:  1 = 11“),  and  a relative  magnetic  permeability  of  unity.  the  incident  wave 


is  polarized  in  the  x direction  so  that  T IT  ,0,0)  and  11  10,11^,11  1. 
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From  (.-.1)  we  then  have  the  following  fundamental  integral  equation  for  the 
determination  of  the  internal  fields  (7eV): 
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lhe  magnetic  fields  are  given  by 
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(2.4b) 


where  the  region  of  integration  V is  the  first  quadrant  (-°°  < x < i®, 
y > 0,  2 > 0). 

A solution  to  (2.5)  is  sought  in  the  following  form,  which  consists 

ot  a superposition  of  a plane  refracted  wave  and  an  unknown  wave  from  the 
★ 

edge  (Aka)  : 


i(k,y+k.z)  ity+ist 

F-X(r)  = Ac  ■’  + — — f(t)dt  (2.5a) 

a 


Aleksandrova  and  Khithnvak  use  the  same  functional  notation  for  the 
weighting  functions  in  t-  and  n-planes  which  is  potent i a l 1 > confusing, 
he  denote  the  weighting  function  in  t-plane  by  t'(t)  to  distinguish  it 
from  ffu). 


Wo  now  write 
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ilf ' = Sj  (p, ip)  + S,(p,<|>) 


where , after  integration  over  V (Appendix  1),  we  obtain 
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So  far  we  are  in  full  agreement  with  Aleksandrova  and  Khithnyak. 
However  it  must  be  pointed  out  that  the  expressions  for  Sj  and  S,,  as  given 
by  ( 2 . Sb ) and  (2.8c),  are  true  only  if  the  following  conditions,  (Al.l) 
and  ( A 1 . 6 ) , are  met  for  ^ real  and  tea. 


tm  k , > 0 
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This  condition  is  not  met  on  a.  However  as  shown  in  Appendix  1\  the 
expression  (2.8cl  tor  S,  is  still  correct  provided  that  the  contour  x is 
deformable,  into  a contour  on  which  this  condition  is  met,  without  crossing 
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5 Frans  format  ion  to  £-  and  n-  planes 

The  following  transformations  are  introduced  which  convert  the 
and  t-plane  integrations  into  contour  integrals  in  and  n-  plane 
respect i ve 1 y , 


ul  = k COS 


(2.  10a) 


v = k sin  £ 


i oh  i 


t = nk  cos  n 


(2.  1 la) 


= nk  sin  n 


(.2.1  lb) 


giving  equations  (Aid))  and  (AklO) 
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The  term  sm  is  missing  from  '.-integral  over  Fj  in  (\K101 


■I-  13 


whore 


t'tn)  = = t\nk  cos  n) 


(2.12c) 


and 


k , = nk  cos  v 


(2.  13:0 


k.  = nk  sin 


(2.13b) 


y = o cos  £ 


: = P sin  f 
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Tito  contours  I'p  and  G are  shown  in  bigs.  3a  and  3b  and  (1  is  shown  in 

(-  IT  / J + i CO  A 

“ I , are  defined 

tt/2  - i'xV 

below  (see  Append i \ II  for  a discussion  about  the  choice  of  the  contour  G) : 

f.  = cosh  ' il  cos  f.  1 ; sin  sinh  — 0 (2.15al 
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n o 

nk  cos  n = t € a when  H G (2.15c) 

o 

l'he  shaded  portions  in  l:ij;4. 3a  and  I ig.4.3h  correspond  to  the  regions  where 
Im  cos  f.  s 0 and  In  cos  (f.  - f ) x 0 respectively. 


•t.J.  I Deformation  of  contours  in  the  .-plane 

Ihe  crucial  step,  in  Aleksandrova  and  Khichnyak's  method,  of  obtaining 
a singular  integral  equation,  involves  the  closing  of  the  contour  at 
infinity  and  deforming  I , onto  C.  In  what  follows  we  will  make  a careful 
examination  of  these  steps. 

4 . 2 . -1 . 1 Reduction  of  integrals  in 

let  us  consider  Sj  as  given  by  (.d.ldal.  According  to  Aleksandrova  and 
khithnyak  the  first  term  in  gives  rise  to  poles  at  cos  •”  = n cos  v and 
sin  . , - n sin  , and  the  second  term  in  Sj  gives  rise  to  poles  at  cos  ^ = 
n cos  , and  sin  . , - n sin  y.  They  seemingly,  conclude  that  the  residues, 
at  sin  A,  - n sin  y,  in  the  tirst  and  the  second  terms  of  Sj  cancel  each 
other,  and  write  down  the  following  expression  (AK11)  for  Sj, 
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whore  the  first  and 
cos  n cos  I,1'  in 

c 1 o ira  t hat  i f A and 


second  terms  in  i-2.il>)  correspond  to  residues  at 
the  first  and  second  terms  of  (,2.12a).  further, 
i|i  are  given  by  [equations  (AK18)  and  (AKlll) 


t hev 


n sin  i|i  *■  si  n >J)  ^ 

n cos  4-'  cos  »|)  (2 . i h 1 
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then  the  first  term  gives  rise  to  a wave  which  is  exactly  equal  to  the 

refracted  wave  A e.\p[ik,y  + if.:.)  and  the  second  term  gives  rise  to  a wave 

which  exactly  cancels  the  incident  wave  f^e.xp  | i k , + ik_^.:  |.  However  there 

are  several  inconsistencies  in  their  statements,  as  described  below. 

i)  With  the  time  dependence  of  the  form  exp [ - itot | the  correct  expression 

for  the  incident  wave,  in  the  fourth  quadrant , is  given  by  exp[-ikpCos  I 

with  .in/2  < J>  < 2 ii  . Since  from  (.2.7)  Imv  '•  0,  the  square  root  expression 
o 

k vl-  n“  cos"  Ip,  in  (,2.1o)  must  be  interpreted  such  that  it  iias  positive 
imaginary  part.  This  in  view  of  C—  - I b)  and  the  tact  that  sin  ■ ■ 

negative  i ni[i  lies  that  k 1 11“  cos  tjh  -k  sin  t.onsequent  ly  , the  second 

term  in  (.2.  In)  takes  the  form 


constant  x exp  [ + ik.>  cos  (.i|>  <1<  ) ) 

which  cannot  cancel  the  incident  wave. 

ii)  When  k,  and  k.  are  given  bv  (.2.12)  |equation  AKS]  the  tirst  term 
,*> 

in  (2 . 5a)  [equation  l AKS)  ) becomes  A exp  [tinkp  cos  (4  40  1 which  i s not  t)ie 

correct  expression  for  the  wave  refracted  at  the  lace  4>  0.  I he  correct 

expression  must  have  a negative  sign  in  the  exponent. 


This  might  lead  us  t o 


I 1') 


suspect  tli.it  the  t Line  dependence  tliat  they  used  is  actual lv  of  the  form 
*■  L u)  t 

e . However  this  will  not  account  for  the  discrepancy  ti),  rey.ard  i ny  the 
extinction  wave.  Also  the  tar  tield  diffracted  wtive  that  results  from  their 
solution  has  a space  variation  of  the  form  exp  [+inkp  + in,  !|/\/nk,  which 
wou'd  he  an  incomint;  cylindrical  wave  if  a time  dependence  of  exp  |*-i..t  | 
is  assumed.  More  importantly  the  type  of  Green's  function  used  in  f d . a ) 
implies  a time  dependence  ot  e l"t.  Hence  we  rule  out  the  possihilitv  of 
the  exp  | M.ot  | t i iik*  dependence  and  conclude  that  their  expression  for  the 
refracted  wave  is  incorrect. 

Since  they  do  not  explicit lv  state  the  values  of  k,.  and  k.,  one  mav  still 

_t)  all 

use  fd.d)  [eipiat  ion  Akd  | to  represent  the  incident  wave  provided  that  k,() 

and  k.,  are  interpreted  as 
ad 

kJ()  - -k  cos  (2.  18a) 

k-0  = -k  sin  IJ.lSbl 

iiil  Matters  will  not  he  resolved  hy  simply  takiiu;  k,  and  k.  to  he 
Icontrtiry  to  equation  AkS| 

k,  = -nk  cos  t|i  (2.l;>a) 

k . = - nk  s i n 4’  id . 1 Phi 

because,  with  and  hence  tji  being  in  the  fourth  qu;idr;int  , k,  will  now  have 
a negative  i magi  nan  part  which  is  in  violation  of  a necessarv  condition 
id. Pa)  for  obtaining  id.  8b).  Nevertheless  k,  and  k.  must  he  given  hv  id.  IP) 
so  that  h lr)  as  ;;  i ven  h\  id.Sa)  will  contain  the  correct  refracted  wave. 


Later  in  section  4 • a we  show  a way  of  getting  around  these  difficulties, 
which  involves  the  decomposition  of  the  v-integral  in  into  two  integrals. 

We  will  find  that  the  result  would  not  only  contain  the  correct  refracted 
wave  hut  also  the  correct  extinction  wave  as  well. 


*•2.1.2  Reduction  of  integrals  in  S, 

Aleksandrova  and  khichnyak  write  down  the  following  expression  for  S, 
[equation  AK15]  without  any  clear  explanation,  except  saying  that  "W'e  consider 
the  singularities  in  S,  in  similar  fashion...", 
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C2.20) 

where  in  the  last  integral,  as  pointed  out  by  them  correctly,  the  order  of 
integration  is  interchanged  in  accordance  with  a corollary  to  the  Poincare 
Bertrand  formula  [34]. 

Their  content  ion, seemingly , is  that  when  P ^ is  closed  at  infinity  and 
F,  is  deformed  into  (1  the  first  double  integral  in  (,2.12b)  gives  rise  to 
two  residues  at  the  poles  n cos  p cos  g and  n sin  n sin  f.  and  the  second 
double  integral  gives  rise  to  a residue  at  the  pole  n sin  p sin  g and  a 
singular  double  integral  which  is  the  second  term  in  (2.20).  The  first  integral 


The  negative  sign  in  front  of  I it  is  missing  in  their  paper,  which  we 
believe  is  a typographical  error. 


Appar- 
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In  (.2.20)  corresponds  to  t ho  residue  of  the  pole  n cos  rj  = cos  £. 
entlv  they  conclude  that  the  residues  due  to  poles  at  n sin  n = sin  \ in 
the  first  and  the  second  terras  of  (,2.12b)  cancel  each  other.  However  a 
careful  analysis,  given  later,  shows  that  this  is  not  true.  Further,  since 
C is  a singular  path  (.cos  £ = n cos  n when  n c and  ('  c C)  one  has  to  give 
due  regard,  to  the  pole  cos  £ = n cos  p,  when  deforming  F,  onto  G.  In 
point  of  fact  one  must  add  a term  corresponding  to  half  a residue,  due  to 
this  pole,  so  that  the  remaining  singular  double  integral  may  be  legiti- 
mately interpreted  as  a Cauchy  principal  value  integral. 

In  the  next  section,  after  obtaining  the  correct  expression  for  S j , 
we  go  on  to  consider  the  above  mentioned  points  in  detail  and  show  that 
the  expression  for  S,  as  given  by  (2.20)  is  incorrect.  The  correct  expres- 
sion for  S 0 will  be  shown  to  contain  additional  terms  one  of  which  corres- 
ponds to  half  a residue,  due  to  the  pole  cos  £ = n cos  n,  and  the  rest  of 
the  terms  arising  due  to  the  non-cancellation  of  the  residues  corresponding 
to  the  poles  sin  = ± n sin  ri  in  the  first  integral  and  the  pole 
sin  t,  = n sin  r)  in  t he  second  integral  in  (2.12b). 
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4.3  couuhc n:.t>  i:\.M,U:vno\  oc  nu  voi u\u  i ntkqui 


Consider  the  volume  integral  (2.Sa)  which  is  reproduced  here  for  coinen- 


E(r') 


i h | f - v ' 


(3.11 


dr  = s | +■  s , 


V |r-r  | 

where  V is  t he  volume  of  t he  dielectric  wedge  and  h(r)  h (r)  >s  given  by 
(,2.5a).  Sj  and  S,  represent  t he  cent  ri  hut  ions  to  the  integral  from  the 
first  and  second  terms  of  I-  in  (2.5a).  Hence 


S,  = A 


i ( k , y ' ' + k . ^ i k | f - f ' | 


(3.  'al 


and 


S , = 


1 1 v + i s c 


f C t Kit 


l , 


ik| f-f ’ | 


dr’ 


(3.2b) 


v U 


Now,  we  wi 11  consider  Sj  and  S,  separately. 


*•3.1  Herivation  of  an  expression  for  S|  containing  the  correct  refracted 
and  the  extinction  waves 

keeping  the  discussion  in  the  previous  section  in  mind  let  us  reaffirm 
that,  with  an  assumed  time  dependence  of  the  form  exp  | - i c t ] , and  with 
EQ(r)  and  l:x(r)  being  given  by  (2.2al  and  (2.5a)  respectively,  the  constants 
k,Q.  k.y,  k,  and  k.  must  satisfy  equations  (2.1S1  and  (2.19)  which  are 


reproduced  here. 


k -> . » = -k  cos 


(3.5a) 


k30  = 'k  sin  *o 


13. 3b ) 


k,  = -Ilk  Cl 


(3.4a) 


k3  = -|,k  si»  <J»  (3.  |b) 

whoro  A anJ  * are  rclateJ  to  4>o  through  (2.17),  and  with  the  implication  that 
V as  shown  in  I-ig.I.i  is  in  the  fourth  quadrant  (_3tr/d  < < dit). 

Now,  we  write  Sj  as 


7 7 ik,y'*ik  " iklr-r'l 

Sj  = A dy ' d;'  e - ''  ° dx' 


r r i k ,v  + i k 

A dv’  I d:'  e 


{IT  i II1  Ik/ly-y'  )“+  ) t 


which  may  be  split  into  the  following  sum. 


S1  = S10  + S11 


(3.  (>a) 


ik ,y ' +ik  . c ' 


s , . , ~ A dy'  dc'  e ' •'  i i II1 1 1 


{ n i *l0  lk/(v  i ' )“♦(.:  ')“  ) 


(3.0b) 


f f i k ,v ' + i k_: ' /- — — — 

Sn  = -A  dy’  j da'  e " {ni  11^ 1 ’ (k.'t.x  y ’)“+(.  = ^ } 

0 (3. be) 


Ill  Appendix  III  t hose  integrals  are  shown  to  he  given  by 


-iko  eos  (**>-({> ,)  , -inkp  cosl£-i|i)  > 

e o + \ o 


k~  in'- 1 ) 


(3.7a) 


i . . . 

' r uev+ 1 k - , 

S..  * iA  \2  [ - — ■ \ : 

1 J ik,-toHk:+ur-k“) 


u*'v  + i v: 
c 

V i k , -ol)  (k  v 


(3.  7b) 


Lot  us  note  the  physical  significance  of  decomposing  into  two  parts  namelv 

Sj^  and  Sjp  The  term  S.^  represents  the  integral  over  the  half  space 
It'  '■  0)  and  Sj j represents  t he  negative  of  the  integral  over  the  quarter 
spa  e ly ' < 0,  t'  > 0).  If  we  were  to  consider  the  problem  of  plane  wave 
incidence  on  a dielectric  filled  half  space  the  integral  equation  for  the 
fields  in  the  dielectric  half  space  would  still  be  the  same  as  given  by 
Id. 3)  with  the  range  of  integration,  now,  being  the  entire  half  space 
It'  '■  0) . Since  we  know  that  the  fields  inside  the  dielectric  are  completely 
given  by  the  single  term  A exp  (ik,v  + ik.tj,  with  k,,  k.,  A,  and  being 
given  by  (5.4)  and  (d.17),  we  expect  that  the  following  equation  must  hold. 

ik  ,v+ik_t  . , . . , 

v e - -•>  = c-iko  COSl4>-4>0) 

'o 

, d d (■  i k , y ' m k . z ' i k | f - f ' | 
k In  - 1 ) I , d o e , . 

+ , \ e i -rr-  d f 

4:i  J I r - r | 

Z>0 

1 A 

-ikp  cosu*)-^,)  k'yr  11 
o In  ‘ 10 
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It  is  a trivial  matter  to  verify  this  if  wo  substitute  for  SJ(1  from  (3.~a) 
and  note  that 


k,y  h.t  = -nke  cos  - y 


l'hus  by  decomposing  S ^ we  have  simply  separated  the  scattered  fields  corres- 
ponding to  the  half  space  problem. 

Let  us  now  consider  S ^ as  given  by  (3.?b).  The  first  integral  on  the 
right  hand  side  may  be  evaluated  by  closing  the  c-cvmtour  with  a large  semi 
circle  in  the  upper  half  plane.  In  the  second  integral  we  may  deform  the 
u)-eontour  into  a branch  cut  contour,  x as  shown  in  Pig. 1.5,  with  due  consid- 
eration to  the  pole  locations.  The  poles  at  - k , - k cos  f and 
1 - o 

2 2. 

to,  = -k  y (.  1 - n“l  + cos“  f lie  below  the  real  axis  and  do  not  contribute 
- o 

to  the  integrals.  Ihe  pole  at  + kV  ( 1 - + cos'  e gives  rise  to  euual 

0»  O *• 

and  opposite  residues  and  hence  does  not  make  any  net  contribution  to  Sjj. 

I lence  S . may  he  s i mp  1 y w r i 1 1 en  as 


sn  - -iA 


l .ON  ♦ W’Z 


V (k  , -u» ) t k . - V ) 
2 .*> 


(5. 10) 


which  after  transforming  into  the  'plane  becomes 


ike  cos  ( l f) 


k“  (n  cos  i)>  +■  cos  in  sin  , +<iu 


(5.111 
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One  »!uy  .is  well  transform  the  u»- intoijra Is  in  (.3.710  into  the  .\-plane  and 
then  make  p rover  contour  detonnations  and  arrive  at  the  same  result.  In  an> 
event  may  finally  be  written  as 


> 

s,  .V  L e ikp  CM«  ♦ \ e ink(  Cl 


k " (,n~-  l ) 


t ko  cos  i - ,*■) 


kJ  ^ 


cos  ) (.n  sin  if  *■  sin  5) 


4.5.J  Reduction  of  the  volume  integral  S,  and  the  derivation  of  the 
singular  integral  equation 

In  sec t i on. I.d . 1 i.e  obtained  an  expression  for  which  contained  the 
correct  form  of  refracted  wave  and  the  extinction  wave.  first  we  noted 
that  k,  and  k.^  must  be  given  by  id. -11  so  that  the  first  term, 

\ exp  [ik,y  ik-c]  . in  id. sal  would  represent  the  correct  refracted  wave. 
Since  k,  as  given  by  (3.  la)  violates  a necessary  condition  [d.i'al  for 
obtaining  (.d.Sbl  we  had  to  go  through  a different  route  which  involved  the 
decomposition  of  S ^ into  two  parts.  This  not  only  resolved  the  above 
difficulty  but  also  resulted  in  the  correct  forms  of  refracted  wave  and 
the  extinction  wave.  However  the  expression  fd.Scl  for  S,  is  still  correct 
as  it  stands  (.Appendix  IV 1 because  the  a- cent  our  is  in  the  upper  half  plane 
of  the  complex  variable  t and  hence  the  requirement  Imt  s 0 is  always  met. 
After  transforming  the  integrations  into  a-  and  q-planos  S,  is  given  by 


(.d.ldb.cl  which  max  be  rewritten  as 


r 


4-d8 


S = S 


:i 


1.5. 13a) 


whore 


s = - ‘ j t'midn  o",k- 
k"  • 


i nkp  sin''  sin  . exp  [ iki  eo  cos  j s i n ' .1 

j in  cos  n-  COS  ~.1  (n'sin'r)  - sin';.) 

* 1 (.3.1 3b ) 


- 1 
k" 


f in)  dp 


exp  I : kf  COS  ) )d 
In  cos  n - cos  f, ) in  sin  n - sin  ) 


1 5 . 1 3c  1 


f (p)  = fir)  = f(nk  cos  n)  (3.1 3 J ) 

However  a careful  consideration  of  the  loci  of  the  poles  cos  n cos  ii 

and  sin  f.  = t n sin  p,  as  n is  varied  on  G (=G  + G I,  shows  that  the 

O 0+  o- 

expression  (d.ddl  for  S , is  incorrect.  In  l;i p.  ).3b,  G and  i.  denote 

respectively  the  parts  of  G that  lie  in  the  upper  and  lower  half  planes  of 

The  contours  G and  G in  l-'ic.l.l  have  similar  mean i nr  with  recard  to 
o + o - 

G in  the  p-plane. 
o 1 

l et  us  first  consider  S,..  When  p £ G the  locus  of  the  pole 

- 1 0*-  1 

sin  ",  n sin  ii  is  chon  hv  I’  + dir  in  1-'  i o .4.3a  while  for  peG  it  is  civon 

+ ‘ 0 - 

by  F + dn.  Since  F + dt  is  outside  the  contour  l'(  the  only  contribution 
to  S,j,  from  the  pole  sin  f,  = n sin  n,  comes  when  peG^.  In  a similar 
manner  we  note  that  t he  pole  sin  ",  n sin  n contributes  to  S , ^ onl\  for 
H e Go  . Hence  the  total  contribution  to  S , ^ from  t ho  term  isin"  n“  sin"  p) 


is  civen  In 


| °M'l  cos(C1-4>)  ]f  Cr  Jn  , f expfikp  cos  z,  + Jf(i  . 

i-  l tcos  i -»  cos  n)cos  7T  J [cos  . i cos 


[cos  . j -ii  cos  ’ 1 ) c o s ' 


13.14) 


where  | = sin  ^ [n  sin  p) 

However  the  pole  cos  •,  = n cos  •;  lies  within  1' 
‘ I is  closed  we  obtain  the  following  result 


. lor  all  *i  £ 1 1 . Hence  when 
i o 


k“(n“-l) 


'v'[  ink.  cos  pi  - ,i  1 | t'[n)d-; 


f f[n)e.\p[ikp  cost'.j-^))] 
k ~ G (cos  ; -II  cos  n)co5 

o>  1 


t'(n)e\p[  ike  cos +40  ]d- 
(cos  -n  cos  n) cos  \ , 


[3. 13) 


1°  consider  S,,  let  us  revert  to  Fig.j.ob.  The  locus  of  the  pole 
sin  \ = n sin  n is  given  by  r [=  T+  + TJ  for  and  for  the  pole 

cos  C n cos  n it  is  given  by  G (=  G+  * G_).  When  1\  is  deformed  onto  i 
°nl>  the  lower  segment  I_  is  crossed.  Hence  the  contribution  due  to  the 
pole  sin  = n sin  n is  given  by 


h 


- exp i i k co >(.C j ]f tn  lr 

k“  l (cos  C.-n  cos  n)cos 

V 1 1 

where  ^ is  again  given  by  (3.14). 

Since  G is  the  locus  of  the  pole  cos  v = n cos  rj  we  cannot  bring  F,  into 
coincidence  with  u without  considering  this  pole.  IV e mav  first  write  s> ^ ^ 
as 


, exp[ikp  cos(C  -4>)]f(n) 

s,,  = — dn 

k“  /„  (cos  ; -n  cos  n)cos 
t>  1 1 

o- 


expfikc  eos(.\-M| 

Clfc, 

(cos  f,  -n  cos  n)  (.sin  t,  -n  sin  q) 

(3.  It 

where  is  a very  small  but  finite  quantity.  To  convert  the  last  term  in 

(3. lb)  into  a "Cauchy  principal  value"  integral  we  must  add  a term  equal 

* 

to  lull:  a residue  due  to  the  pole  cos  " = n cos  q resulting  in 


f (n)dq 


Ibis  was  related  to  Aleksandrova  and  Khizhnyak  in  a private  communication 
lo  date  there  has  been  no  response  from  them. 
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Vr  j-  exp  [ Lk.  cos  f.  j-4*)  ] f (n  )Jp 
. 2 J (cos  ; . -n  cos  r)  Jco  . 


exp  I i k.'  cos  (C  f (n ) dn 

sin  £ ->  ( s in  ",  , - n sin  n I 


, . , - nl  f (n > dn 

explikp  cos  l |d  . = . ^ . , , ■, — 

■ 1 J (.COS  -II  COS  nits  in 


- n sin  n) 


(5.  17:0 


= cos  (n  cos  n) 


(5.17b) 


whore,  in  the  last  term,  we  have  changed  the  order  of  integration  which 
is  permissible  according  to  a corollary  to  the  Poincare-Bertrand  formula, 
lhe  p- integral,  in  the  last  term  of  (3. It's),  may  now  be  correctly  inter- 
pretted  as  a "Cauchy  principal  value"  integral.  Combining  (3.15)  and  (3.17) 
we  arrive  at  the  following  expression  for  S,. 


exp[  ink.'  cos (tt -4> ) ] f (n)dn 


k~ (n“- 1 ) q 


f (n)dn 


x]i[ikp  cost*.  'M  Ids  tcos  t>n  cos  n)(sinh-n  sin  n) 


exp[ikp  cos ( ",  j - ) 1 f (n ) dn  ...  r exp [ ikp  cos (£j+  i 1 ] f (P) dp 


(cos  c^-n  cos  n)cos  i j (cos  gj-n  cos  p)cos  r, ^ 

G 

o- 


,'xp[  ikp  cos(r ; -if ) | f tn)dn 


exp[ikp  cos  (7  ,-ij> ) | f (n) 


,2  j (cos  £.-n  cos  nlcos  C.  ,2  J sin  :;,(sin  7 -n  snip  1 

k G 1 K (1 

o-  o 


(3. 18) 


winch  has  four  additional  terms  as  compared  with  the  expression  ( \ K 1 S l 
obtained  hv  Aleksandrova  and  kh  i dunak  for  s . if 

substitute  for  Sj  and  S,,  which  determine  the  volume  integral  in  i'.o) 
fa. Id)  anil  (a. IS)  we  obtain  the  following  integral  equation  for  flu). 


e\p|ikp  cost.',  i|l)|d.. 


\ 

i|»)  fs in 


Ml  M II  , I 1 (cos 

l! 

o 


:>  i 


f exp  | i ke  cos  ( v ) | f (n  )dtl 
(COS  | -1!  COS  n)COS  «' 


.':i  i 


exp  [ i ko  COS  f * ;■  I I f ri  ■ dt 

(cos  n cos  q)  cos 

; 1 

o 


r exp  | iko  cosih  -iji)  | f(n)dn 
1 j i cos  ■ , n cos  i) ) cos 


t e\p | ikp  cos i r. , ;>)  | f(n)dn 
• n i I 0 . 

s i n . . , (si  n . ,-n  sin  q ) 


where  ^ and  are  defined  b>  (a. 11)  and  (a. 17b).  If  the  last  four  t 
were  not  present  in  (a.  Id)  the  resulting  equation  could  be  solved  exac 
by  using  standard  techniques  [alpah|,  as  has  been  done  by  Aleksandrova 
kin  hnyak.  We  see  no  way  of  avoiding,  the  presence  of  these  additional 
.uni  hence  conclude  that  Aleksandrova  and  Kh i; hnyak’ s method  fails  to  le 


, from 

f (M Id" 
i cos n)i 


Id. Id) 

e rms 
tlx. 
and 
t e rms 
ad 


to  a tractable  integral  euqat ion  for  the  weighting  function  f(q). 


i-  n 


-t.  i I'OKMiii.vi'iON  usim;  a Manmr.n  contour 

Aleksandrova  and  Khichnyak  seek  the  solution  for  the  diffracted  fields 
in  the  form  of  an  integral  over  a branch  cut  contour  a.  It  was  felt  that 
such  a representation  of  the  diffracted  fields  might  not  he  complete.  \n 
integral  over  the  real  axis  seemed  to  he  more  appropriate.  With  this  idea 
in  mind  we  could  seek  a solution  in  the  following  form 


i tk  ,y  < k ) 


r x 1 1- ) \ e 


ttvHs: 

e 


f It  hit 


pi  .la) 


where,  as  in  previous  sections 


k“  + k“  = n“  k" 


(4.2a 


t“  + s“  - n “ k’ 


1 4.2b) 


k,  = -Ilk  COS  l/l 


14.5a) 


k . = - ii  k s i n tit 

o 


(4.  oh) 


cos  iji  (cos  ^ 1 / n 
o 


(4.4a) 


A 


21  sin  41 

o o 

non.  * sin  | 


( 1.  Ih) 
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It  is  assumed  that  the  incident  wave  is  in  the  fourth  quadrant  as  given  by 
(2.2)  while  k and  k.()  are  defined  by  (.3.31,  and  as  shown  in  I'iyq.n  t !te 
contour  ok  integration  t'  ( * •!'})  in  t lie  t-plane  is  chosen  to  lie  just 

above  the  real  ix:  . i\  ^ 01.  In  the  limit  we  let  A k 0 with  the  under- 

A 

standing  that  am  poles  of  t"  (.  1 1 that  lie  on  the  real  axis  of  t are  indented 
with  clockwise  sev.ii -c  i rc  les . Ke  split  the  volume  integral  appearing  in 


(J.31  into 

the 

sum 

S j ♦ S where  S ^ 

. S,  are  given  by  i.3..'). 

In  the  sum 

S S ♦ 

1 10 

11 

t ho 

express  lens  for 

S and  S ^ are  given  by 

(3. "1  which 

are  reproduced  here 


10 


111 

k“  (n“- 1 1 


j:  -ikp  C0S(M>01  + N o-inkp  cost;  ,)  k t'.-'l 


Sl, 


iojv+  \ k~z 

j ' > 7 d 

( k , - 0 ) - k k . 1 


iuv  t iv: 
o 

v . . ' ( k . \ i 


do)  r e.e 


The  first  integral  in  Sj|  may  be  evaluated  by  closing  the  .'-contour  with  a 
large  semi  circle  in  the  upper  half  plane  and  in  the  second  integral  the 
path  of  the  integral  may  be  shifted  slightly  above  the  real  axis  with  the 
following  result. 


i v k -k:  y + i k . 

-2ttA  e 

Ik  , v^k-  -k“  ) /k“-k“ 


- i.\ 


In  writing  (.1.3)  we  note  that  the  pole  u> 
hence  does  not  contribute  any  residue. 

ero  such  that  the  above  deformation  in 
without  crossing  the  poles  at  \ k“ 

except  for  the  special  case  when  $ 3 

cons i derat  ion  here . 


i wv  + i v ~ 

c d.  (4. 

V ( k , -tO  ) ( k . - V ) 

k,  is  below  the  real  axis  and 
We  choose  A sufficient 1>  close  to 
the  second  integral  is  possible 
kt  . fit  is  is  always  po  sihle 
which  we  may  exclude  from 
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Proceeding  in  .1  manner  similar  to  that  in  Appendix  1 we  obtain  t he 
following  expression  for  S , 


S , - 1 


■r,  uuv  + is; 
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L 
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(-l.'dal 


Im(s  + v)  > 0 


l l.'dbl 


which  are  satisfied  for  any  arbitrarily  small  A N 0. 

Hie  first  w- integral  in  f f . S 3 may  he  evaluated  by  closing  the  contour 
with  a large  semicircle  and  the  second  integral  may  be  converted  into  a 
"fanchy  principle  value  integral”  by  properly  accounting  for  a hall  residue 
corresponding  to  the  pole  at  to  = t.  When  this  is  done  S,  is  given  by 
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where  the  square  root  expression  1 k"  s'  is  to  be  interpreted  such  that 


1 m>  h ’ 


0 , t £ C. 


I 1 . 1 11 
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If  wo  substitute  for  Sj  and  S,  in  U\3)  we  obtain  the  followim;  integral 


oiiuation  for  t'(t). 


Hoy  * i v: 
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1 1 . l : i 


Because  of  the  presence  of  the  last  three  terms,  which  are  not  self-cancell 
*ni'’  '*  ' s not  possible  to  solve  f 1 - 1 - ) (or  the  weight  in;;  function  f(tl  1" 
UsinR  tho  standard  technique  that  was  employed  by  Aleksandrova  and  khizhnyak 
to  solve  oqn  IAK13).  Tluis . choosing  the  contour  of  integration  alone,  the 
real  axis  does  not  seem  to  avoid  the  difficulties  under  Ivins  the  method  of 
Aleksandrova  and  khizhmak. 
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4.5  CONCLUSIONS 

In  this  report  wo  have  made  a detailed  examination  ot'  Aleksandrova  and 
khichnyak's  method  led]  of  obtaining  a solution  to  the  problem  of  electro- 
magnetic plane  wave  diffraction  by  a rectangular  dielectric  wedge.  P> 
carrying  out  the  spatial  integrals,  involved  in  the  problem  in  full  detail 
and  by  a careful  look  at  the  subsequent  contour  deformat i ms,  we  Have 
shown  that  Aleksandrova  and  khichnyak,  besides  mixing  up  the  sign  convention 
in  the  time  variation,  omit  some  terms  which,  if  properl>  accounted  for, 
would  lead  to  an  integral  equation  that  is  not  only  different  from  what 
they  have  finally  solved,  but  is  also  not  amenable  to  solution  by  presently 
known  techniques.  We  have  attempted  a formulation  using  a modified  contour 
of  integration  in  the  t-plane  but  failed  to  obtain  an  integral  oquat:  n that 
can  be  solved.  We  conclude  that  not  only  is  the  solution  given  by 
Aleksandrova  and  khirlmvak  incorrect  but  also  that  their  method  as  presented 
is  not  capable  of  leading  to  a tractable  integral  equation  for  the  unknown 
weighting  function.  Thus  the  problem  of  wave  diffraction  by  a dielectric 
wedge  remains,  is  yet,  unsolved,  even  for  the  special  case  of  a rectangular 
wedge,  hut  for  the  solution  of  Uadlow  |10)  which  has  been  questioned  by 
others  [20,21]. 

The  solution,  of  Aleksandrova  and  khi  hnyak,  to  the  problem  of 
diffraction  by  an  arbitrary  angled  dielectric  wedge  [55],  is  also  open  to 
debate  since  the  results  are  said  to  he  based  on  a method  identical  to  the 
one  used  in  their  earlier  paper  [52]  which  we  have  discussed  in  this  report. 
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for  S,  subject  to  the  following  condition: 

Im  t > 0 (Al.ba) 

Im(.s  H')  > 0 when  t e a,  and  a e real  axis  (Al.bb) 

However  since  s takes  on  different  signs  on  two  sides  of  the  branch 

cut  tlie  inequality  (A1.6b)  is  violated  on  part  of  the  contour  a . To 

avoid  this  difficult)'  we  do  the  following: 

We  deform  into  another  contour  a'  where,  as  shown  in  Fig.B  1.1. 

the  latter  coincides  with  « on  the  right  side  of  the  branch  cut  and  on 

the  left  side  it  is  defined  bv  the  equation  t't"  = j n'k"  where  t =t'  + it” 
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and  k = k (1  + id).  We  may  now  carry  out  the  spatial  integration  since 
the  conditions  (.Al.ba.b)  are  satisfied  for  t r a.'.  After  carrying  out 
the  spatial  integration  the  contour  a'  may  be  replaced  by  a provided 
the  weighting  function  fit)  meets  certain  requirements  as  discussed 
in  Appendix  IV. 
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Aleksandrova  and  Khicnyak  do  not  explicitly  state  what  G is  except 
indicating  that  G is  ( + |^)  which  only  specifies  the  end  points 

of  G.  In  their  later  paper  [35]  on  diffraction  by  an  arbitrary  dielec- 
tric wedge  the'  explicitly  define  G in  the  rj-plane  (not  the  contour  G 
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the  exact  shape  yet  to  be  determined.  We  first  consider  G to  be 
given  by  (2.15a)  which  leads  to  vertical  branch  cuts  in  the  t-plane  as 
shown  in  Hi  g. 4.2.  This  necessitates  some  manipulations  (Appendix  IV) 
in  obtaining  the  expression  (2. Sc)  for  S,  so  that  the  necessary  conver- 
gence criteria  are  met  and  also  places  some  restrictions  on  the  singu- 
larities of  f(t).  If  G is  taken  to  consist  of  straight  line  segments 
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and  pass  through  the  points  I-tt/2  + i , C_Tr/  — + 10).  C71  / - + 10)  and 
/ 2 - l“>) , then  the  corresponding  contour  a in  t-plane  would  be  such 
that  lin  s = 0 for  tea  and  I in  s > 0 in  the  entire  Riemann  sheet. 
This  necessitates  a modification  of  the  shape  of  branch  cuts  in  the  t- 
plane.  liven  though  such  a choice  of  branch  cuts  seems  more  natural  for 
our  problem  our  attempts  17  using  such  a modified  contour  also  failed  to 
give  a desired  singular  integral  equation  for  the  weighting  function 
fit). 
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The  exponential  integrals  may  now  he  evaluated  to  give 
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subject  to  the  following  conditions 


Im(k.i-v)  s 0 
III!  k , < <) 


(A3. 15a 
(A3. I Sb 


which  are  clearly  satisfied. 


liquation  (AA.ll)  is  the  same  as  (3.7bl. 


AI’PIINPIX  IV 


V. VALUATION  0i:  I'll!  I XiT.GRAI.  S, 


In  this  appendix  wo  derive  the  expression  (2.8c)  for  S,  which  is 
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where  V e.  (-00  <x  <»,  y >0,  z N0)  and  a is  shown  in  Tig.  B4.1. 

We  note  that,  with  the  choice  of  a ns  shown  in  fig.  Rd.l,  s is 
complex  and  has  opposite  signs  on  the  two  sides  of  the  branch  cut.  Ilenct 
the  exponential  term  expjity'  + ist'l  will  diverge,  as  vo°,  on  one 
side  of  the  branch  cut  if  it  is  convergent  on  the  other  side  and  vice 
versa.  To  overcome  this  difficulty  we  replace  a by  a'  on  which 
1m  s _>  .).  This  is  permissible  if  the  following  conditions  are  satisfied, 
i)  f(t)  has  no  singularities  between  n and  o' 
ii)  f(t)  takes  a finite  limit  as  t v 


The  second  condition  is  needed  to  ensure  that  the  t: -integral  over  the 
segment  AA1  is  zero.  1'hus  the  final  solution  for  f(t)  must  be 
checked  against  these  conditions.  It  may  be  worth  noting  here  that  even 
though  we  believe  and  show  in  this  report  that  the  function  f(t1  = f(til 
given  in  Aleksandrova  and  khi  .-hnyak's  paper  is  obtained  by  solving  a 
wrong  integral  equation,  it  does  indeed  meet  the  above  two  requirements 
keeping  the  previous  discussion  in  mind  we  deform  the  contour  a hi 
(Hi.  11  into  t and  after  changing  the  order  of  integration,  write  S,  a- 
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where  we  have  made  use  of  the  Ilankel  function  representation  (2.o). 

Now  we  may  carry  out  the  y'  and  z'  integrals  by  splitting  the  range 
of  y'  into  y'<  y and  y'  >y  and  correspondingly  the  range  of  z'  into 
z <z  and  z'  >z  and  obtain 
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The  conditions  (A4.4a,b)  are  obviously  satisfied  now.  IVe  note  that  the 
terms  in  the  integrand  that  are  odd  functions  of  to  , give  zero  after 
integration  over  the  real  axis  and  hence  S-  may  be  written  as 
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Now  by  a reasoning  similar  to  that  used  before  we  may  now  deform  a' 
into  a and  obtain  (2.8c). 


